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The  objective  of  this  program  is  the  exploratory  development  of  doubly  rotated  cuts  of  quartz  possessing 
superior  Surface  Acoustic  Wave  (SAW)  properties  for  applications  involving  environmentally  hardened 
devices.  The  key  properties  examined  and  optimized  both  theoretically  and  experimentally  are:  first,  second 
and  third  order  Temperature  Coefficients  of  Delay  (TCD),  piezoelectric  coupling  factor,  power  flow  angle, 
Bulk  Acoustic  Wave  (BAW)  inverse  velocity  surfaces,  degeneracies  and  leaky  waves,  and  sensitivities  of  the 
above  quantities  to  misorientations  and  manufacturing  tolerances. _ >_ 
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^The  program  consists  ot  two  major  task  areas  comprising  an  interactive  numerical /experimental 
approach.  Task  I  involves  the  numerical  computation  of  the  key  SAW  properties  for  doubly  rotated  quartz 
substrates  for  the  purpose  of  locating  promising  angular  ranges  with  properties  superior  to  the  singly 
rotated  cuts  now  in  existence.  More  detailed  calculations  follow  to  refine  the  angular  coordinates  in  order  to 
specify  cuts  for  experimental  verification  inTask  II.  In  Task  II,  sets  of  substrates  with  promising  orientations 
identified  in  Task  I  will  be  prepared  and  SAW  device  patterns  will  be  fabricated  for  evaluation  of  the  key  SAW 
properties.  The  experimental  results  of  this  task  will  be  correlated  with  the  theoretical  predictions  and  an 
iterative  process  develops  for  refinement  of  both  theoretical  and  experimental  parameters.  As  the  program 
proceeds,  working  SAW  device  models  will  be  delivered  as  a  demonstration  of  progress  and  an  indication  of 
the  future  potential  of  the  doubly  rotated  cuts  JJepending  upon  the  progress  made  and  time  and  budget 
limitations,  additional  properties  in  the  area  of  nonlinear  elasticity  will  be  investigated.  This  report  contains 
the  result  of  Task  I  and  part  of  Task  II.  \ 
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SECTION  1 
INTRODUCTION 


1.  PROGRAM  OBJECTIVE 

The  objective  of  this  program  is  the  exploratory  development  of  doubly  rotated  cuts  of  quartz 
possessing  superior  Surface  Acoustic  Wave  (SAW)  properties  for  applications  involving  environmentally 
hardened  devices.  The  key  properties  examined  and  optimized  both  theoretically  and  experimentally  are: 
first,  second  and  third  order  Temperature  Coefficients  of  Delay  (TCD),  piezoelectric  coupling  factor,  power 
flow  angle,  Bulk  Acoustic  Wave  (BAW)  inverse  velocity  surfaces,  leaky  waves,  and  sensitivities  of  the  above 
quantities  to  misorientations  and  manufacturing  tolerances. 

2.  PROGRAM  SCOPE 

The  program  consists  of  two  major  task  areas  comprising  an  interactive  numerical /experimental 
approach.  Task  I  involves  the  numerical  computation  of  the  key  SAW  properties  for  doubly  rotated  quartz 
substrates  for  the  purpose  of  locating  promising  angular  ranges  with  properties  superior  to  the  singly 
rotated  cuts  now  in  existence.  More  detailed  calculations  follow  to  refine  the  angular  coordinates  in  order  to 
specify  cuts  for  experimental  verification  in  Task  II.  In  Task  II,  sets  of  substrates  with  promising  orientations 
identified  in  Task  I  are  prepared  and  SAW  device  patterns  fabricated  for  evaluation  of  the  key  SAW 
properties.  The  experimental  results  of  this  task  are  correlated  with  the  theoretical  predictions  and  an 
iterative  process  develops  for  refinement  of  both  theoretical  and  experimental  parameters.  As  the  program 
proceeds,  working  SAW  device  models  will  be  delivered  as  a  demonstration  of  progress  and  an  indication  of 
the  future  potential  of  the  doubly  rotated  cuts.  Depending  upon  the  progress  made  and  time  and  budget 
limitations,  additional  properties  in  the  area  of  nonlinear  elasticity  will  be  investigated. 

3.  TECHNICAL  APPROACH  SUMMARY 

To  accurately  characterize  the  properties  of  doubly  rotated  quartz,  three  basic  capabilities  are 
essential: 

a.  Theoretical  approach  and  associated  computer  software  which  will  accurately  and  quickly 
locate  promising  zero  TCD  cut  angles  and  characterize  the  other  key  SAW  parameters 

b.  Source  of  rotated  quartz  substrates  of  superior  quality  which  can  be  quickly  fabricated  and 
the  angular  orientation  which  can  be  determined  with  a  high  degree  of  precision 

c.  Required  fabrication  facilities  and  measurement  tools  to  accurately  determine  the  key  SAW 
device  parameters. 
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In  the  theoretical  area,  this  program  has  charaterized  two  basic  theoretical  approaches  tor  the 
identification  of  zero  TCDs  on  rotated  cuts  of  quartz.  For  this  study,  two  computer  programs  available  at 
Motorola  are  used.  The  first  program  calculates  the  first,  second  and  third  order  TCOs  of  rotated  cuts  using  a 
finite  difference  method.  This  technique  is  simple,  well  established,  and  has  been  used  for  analytically 
determining  the  temperature  coefficient  curves  for  singly  and  doubly  rotated  cuts  of  quartz.  To  more 
accurately  refine  the  temperature  coefficient  properties,  a  second  program  which  encompasses  lattice 
skewing  effects  is  used.  This  more  complete  theoretical  approach  is  based  on  the  work  of  Sinha  and 
Tiersten1;  its  utility  has  been  verified. 

The  final  theoretical  work  is  the  characterization  of  the  other  key  parameters  with  standard  SAW 
programs  used  routinely  for  material  characterization  and  device  development. 

Accurately  oriented  quartz  bars,  supplied  by  Motorola,  Carlisle,  are  cut,  lapped  and  finely  polished 
at  Motorola.  A  mechanical  polishing  procedure  is  used.  During  this  program,  several  substrates  from  a 
single  bar  with  incremental  angular  deviations  about  a  promising  angular  position  are  fabricated.  By  careful 
organization  of  the  angle  selection  and  cut  procedures,  a  substantial  savings  in  time  and  money  is  achieved. 

The  angular  orientation  of  the  doubly  rotated  substrates  are  defined  to  an  accuracy  of  within  ±5 
minutes  using  X-ray  diffractometry.  Equipment  used  includes  Laue  pattern  X-ray  equipment,  X-ray 
diffractometers,  and  precision  wafer  cut  and  polishing  equipment. 

A  complete  SAW  test  area  and  optical  laboratory  form  the  basis  for  evaluating  the  key  SAW 
parameters  of  the  doubly  rotated  quartz  delay  lines,  oscillators  and  resonators.  The  equipment  will  be  set  up 
for  rapid  display  and  recording  of  the  SAW  parameters. 


1  “On  The  Temperature  Dependence  of  the  Velocity  of  Surface  Waves  in  Quartz",  B.K.  Sinha  and  H.F.  Tiersten 
Proceedings  of  the  32nd  Annual  Symposium  on  Frequency  Control,  1978,  pp.  150  -  153. 


SECTION  II 

TECHNICAL  DISCUSSION 


The  following  section  discusses  the  technical  approaches  used  in  this  program.  A  comparison  of 
numerical  approaches  is  made  and  the  numerical  approach  used  is  described  in  paragraph  1  (Task  I). 
Experimental  techniques  (Task  II)  will  be  outlined  in  paragraph  2. 

t.  INTRODUCTION 

Quartz  is  the  most  commonly  used  substrate  for  fabricating  Surface  Acoustic  Wave  (SAW)  devices. 
In  SAW  narrowband  filter,  oscillator,  and  resonator  applications,  the  temperature  stability  of  the  device  is 
an  important  design  parameter.  Currently,  almost  all  SAW  devices  fabricated  on  quartz  use  the  ST-Cut, 
which  exhibits  a  parabolic  frequency  dependence  in  temperature.  For  many  applications,  the  temperature 
dependence  of  devices  fabricated  on  ST  quartz  is  too  large.  Thus  it  is  desirable  to  find  crystal  cuts  with 
superior  temperature  performance.  Of  course,  many  other  design  parameters  must  be  considered  when 
choosing  a  crystal  cut.  Some  of  the  more  important  ones  are  the  piezoelectric  coupling  coefficient,  acoustic 
losses,  dependence  of  device  performance  on  cut  misorientation,  excitation  of  bulk  modes,  and  beam 
steering  angle.  These  parameters  are  all  easily  determined  for  a  given  cut. 

The  objectives  of  this  program  are  to  find  crystal  cuts  which  exhibit  zero  temperature  coefficients  of 
delay  so  that  there  will  be  no  frequency-temperature  dependence  observed  in  temperature  stable 
oscillators,  resonators  and  filters.  We  have  used  computer  models  to  investigate  the  temperature 
dependence  of  different  cuts  of  crystal  for  SAW  devices. 

Defining  r  as  the  delay  time  for  an  acoustic  wave  to  propagate  between  two  points  on  the  surface  of 
the  crystal,  we  wish  to  find  orientations  for  which  r  is  constant  in  temperature,  or  more  formally  (1/r) 

(dr/dT)  |  1  _  M  tlcftretA  ('  =  0. 

If  F  is  the  frequency  of  a  SAW  resonator,  we  have  (1/F)  (dF/dT)  =  -(1/r)  (dr/dT).  Of  course,  we 
would  also  like  the  higher  order  derivatives  to  be  as  close  to  zero  as  possible,  or  (1  /2r)  (d2r/dT2)|  T . 2<  *,.*«  c 
=  0,  etc. 

Letting  9  be  the  length  between  two  points,  (r)  is  simply  given  by  r  =  (G/Vs),  (1/?)(dr/dT)  = 
(1/?)(d9/dT)  -  (1/V,)(dVs/dT). 

We  have  computer  programs  for  calculating  r  as  a  function  of  the  stress,  the  dielectric  and 
piezoelectric  constants  of  a  substrate  material.  Furthermore,  we  have  at  our  disposal  the  temperature 
variation  of  those  constants  for  quartz  allowing  one  to  calculate  (1  /  V,(dV,/dT)  by  a  finite  difference  method 
(note  that  higher  order  terms  can  be  calculated  the  same  way).  If  the  crystal  expansion  coefficients  are 
known,  it  is  then  a  simple  matter  to  calculate  (1/r)(dr/dT)  for  any  particular  orientation.  Perturbation 
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programs  developed  by  Sinha  and  Tiersten  to  calculate  (1/r)(di7dT)  are  also  available.  Combined  with  a 
search  method,  we  are  able  to  find  cuts  for  which  first  order  temperature  coefficient  of  frequency  vanishes. 


a.  Calculation  of  Temperature  Coefficients 

It  has  been  shown  that  determining  the  temperature  dependence  of  r  (time  delay)  is  equivalent 
to  determining  the  temperature  dependence  of  F  (frequency)  via  the  relation  F  «  i/T. 

Since  our  experimental  data  is  derived  from  frequency  measurements,  we  determine  the 
frequency  characteristics  of  the  devices.  The  relation  between  the  temperature  coefficient  of  frequency 
(TCF)  and  temperature  coefficient  of  delay  (TCD)  are  related  as  follows  (See  Appendix  E): 

^  =  -«/"  (1) 

=  -a  t2)  +  («r'")2  (2) 

=  -a  '"  +  2a';1  a'? -(a ‘I’)3  (3) 

where  ai is  the  ith  order  TCF,  a  T is  the  ith  order  TCD. 

Using  the  relations  above,  one  can  always  relate  one  set  of  the  temperature  coefficients  to  the 

other. 


The  various  procedures  outlined  in  the  following  sections  will  yield  the  quantities  a,1",  i  =  1 , 2, 

3,  with  V,(T)  S  V,(T„)(1  +  a.'1'  dT  +  av(2’  dT2  +  av<3’  dT3)  and  dT  =  T  -  T„. 

The  problem  of  interest  is  not  only,  however,  in  finding  the  a,(i|’s  but  in  finding  the  delay  time  r 
and  the  frequency  dependence  F  of  a  device.  The  frequency  dependence  F(T)  =  F(T0)(1  +  aF'"  dT  +  aF(2’  dT2  + 
a i121  dT3)  is  a  function  of  not  only  Vs  but  of  2(T0)  =  2(T0)(1  +  a,"’  dT  +  a/21  dT2  +  a/3’  dT3),  the  spacing  between 
reflectors  in  a  resonator  (or  the  length  of  the  delay  line  in  an  oscillator) 

From  the  above  relations  and  (1  +  X)-1  s  1  -  X  +  X2  -  X3  it  X  «1,  we  see  that: 

F  =  V,  2  S£  (  V,  (1  +  av"’  dT  +  av(2)  dT2  +  av13’  dT3)}/2.  (1  +  a,*"  dT  +  a,'2'  dT2  +  a,'3'  dT3)  (4) 

sF(To)  { 1  +  av"' dT  +  av,2>  dT2  +  ay,w  dT1  (5) 

-  a,"’  dT  -  a,"*  av'"  dT2  -  a,"'  av'2'  dT3 

-  a,12'  dT2  -  a,'2'  a>'"  dT3  -  a,'3’  dT3 

+  (a,"1)2  dT2  +  (a,"')2  av'"  dT3 

+  2a,'"  a,'2'  dT3  -  (a,"'  dT)3}  (6) 


or 


arm 

=  Qx 

—  Ot  1 

(7) 

arU) 

=  Qfv 

U’1  —  a,"’  ai 

"  +  (a,' 

y  -  a,'2' 

(8) 

Ot\ 

= 

|M  III  1 

—  a  i  av 

21  +  a," 

1  (a,'")2 

(9) 

,"r 

+  2a,"'  a^'1 

t  1  l 

—  Ot 

a,’2’  -  a,1” 

giving  the  frequency  dependence  directly.  The  calculation  of  the  temperature  coefficients  of  velocity,  is 
achieved  by  calculating  the  phase  velocity  with  the  Finite  Difference  Technique  or  other  techniques 
discussed  below  for  a  variety  of  temperatures.  As  the  velocity  is  a  function  of  temperature,  a  linear 
regression  program  is  used  in  the  finite  difference  technique  to  curve  fit  the  data  to  a  third  order 
polynominal.  The  constants  a,"1,  a,121  and  a,13’  are  thus  obtained  by  optimum  curve  fitting  of  the  data  points 
to  the  polynominal.  The  temperature  coefficients  of  length,  a,"1,  are  found  in  standard  references.  It  should  be 
noted  that  the  <*i’s  also  depend  on  the  direction  of  propagation,  and  must  be  calculated  for  each  direction  of 
propagation  considered  by  a  simple  geometrical  transformation  analogous  to  that  used  to  rotate  all  of  the 
other  physical  constants. 

b.  Rayleigh  Wave 

We  will  briefly  discuss  the  Rayleigh  wave  solutions  in  their  general  form  in  this  section.  The 
coordinate  system  is  defined  with  the  Z  or  1  axis  being  the  direction  of  propagation  and  the  -Y  or  3  axis 
normal  to  the  crystal  surface.  V,  will  denote  the  Rayleigh  wave  velocity,  pR = «>/ V„  the  wave  number,  u*  (i  =  1 , 
2, 3),  the  particle  displacements  along  the  1, 2,  or  3  axis,  u4  =  <t>,  the  electric  potential.  Also,  T*  denotes  the 
stress  tensor. 

S„  =  (1/2)(du,/dx,  +  dUj/dx,),  the  strain  tensor,  c,Jki,  the  elastic  tensor,  p  the  density  of  the 
substrate,  D,  the  electric  displacement,  E,  the  electric  field, the  dielectric  constant,  e,Jk  the  piezoelectric 
constant.  We  also  use  a  for  the  difference,  where,  for  example, 

A p  =  p(T)  -  p(T0)  =  p(T0)(o  p("dT  +  a  p<2'dT2  +  a  p131  dT3),  and  let  P  =  power/unit  width  in  the 
x  direction. 

We  assume  relations  such  as 

T  =  — e  •  E  +  cE:S  or  S  =  e-E  +  SE:T 
D  =  <SE  +  e:S  D  =  «TE  +  d:T 
p32  Uj/dt2  -  ciiki  32uk/ax,  dXi  -  ekij  32<*>/0Xi  dxk  =  0 
e,u  32uk/3Xi  3X|  -  e.it  32  4>ldXi  dxk  =  0 


00) 

(11) 
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The  general  solution  for  a  wave  traveling  on  the  surface  can  be  written 


u,= [  S  Cm  a,""’  exp(-(i<u/V.)  2  br  Xi)  ]  expM)  for  j  =  1  to  4 

m  1=1 

For  the  Rayleigh  wave,  this  reduces  to  U,  =  J  cm  exp  (-i/3Rb'm,y)  expflfrz  -  ia>t) 


(12) 

(13) 


where  the  coefficient  cm,  ar’(weighting  factors),  Mwave  number),  and  b,,n'(complex  decay  constants)  are 
to  be  calculated  by  the  standard  iterative  procedure  on  a  computer. 

The  fields  are  calculated  as  E  =  Vu4  with  E,  =  0,  Ey  =  (— i/3Rb,m>)U4  and  E*  =  (i0R)u4. 

c.  Methods  for  Calculating  the  Temperature  Dependence  of  the  Rayleigh  Velocity 

In  the  following  sections,  four  methods  for  calculating  the  temperature  characteristics  of  the 
Rayleigh  wave  are  discussed. 


(1)  Finite  Difference  Technique  of  Calculating  Temperature  Dependence  of  the  Rayleigh 

Velocity 


Before  TCD'"  =  (1/r)(dr/dT)  can  be  calculated,  one  must  first  calculate  the  dependence 
of  the  Rayleigh  wave  velocity  on  temperature.  The  most  straightforward  method  for  doing  this  is  the  finite 
difference  method.  The  Rayleigh  wave  velocities  are  calculated  for  different  temperatures,  yielding  the 

values  V,(T ),  i  =  1 , 2 . n.  This  is  done  by  first  calculating  the  fundamental  constants  at  the  temperature  T<  of 

interest.  The  fundamental  constants  are  then  rotated  into  the  coordinate  system  of  interest.  An  iterative 
procedure  is  used  to  calculate  a  velocity  V,  for  which  Christoffel's  equation  and  the  boundary  conditions 
are  satisfied  simultaneously  (see  Appendix  D).  Simple  finite  difference  techniques  can  be  used  to  calculate 
(1  /V,)(dV,/dT),  (1  /2V,)(d2V,/dT:),  etc.  For  example,  after  calculating  VS(T,)  for  T,  =  To.  T2  =  T„  +  aT,  T,  =  T„  - 
aT,  we  can  use  (1/V,)(dV,/dT)  at  T„  *  (1/Vo(To))(V.(T0  +  AT)  -  V.(T0-  aT))/2a.  Alternately,  standard  linear 
regression  of  polynomials  may  be  employed  to  yield  those  coefficients.  Another  approach  consists  of 
calculating  directly  the  frequency-temperature  characteristics  of  the  orientation  for  several  widely  differing 
temperatures.  A  measure  of  the  temperature  stability  is  then  used.  Letting  F(T)  be  the  frequency  of  the 
device,  the  measure  of  deviation  is  calculated  as 


RMS  frequency  deviation  = 


[{X(F(T.)-F(To))2)/n  F  (T0)2]'  2  with 

i  =  1 

F(T)  =  V.  (T)/2q(1  +  «l'"(AT)  +  oti_(2|(AT)2  +  aL,3'(AT)J), 


(14) 


although  on"',  ah'2’,  and  an13'  are  calculated  for  purposes  of  comparison. 
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In  much  of  the  earlier  work,  the  SAW  velocities  were  calculated  at  three  temperature 
points,  T i  =  -50°C,  T 2 = 25°C,  and  T3 = 100°C,  to  save  on  computer  time.  The  following  formulas  were  used  for 
the  calculations  of  first  and  second  order  temperature  coefficients  of  velocity. 

=  IV(Ti)  -  V(T ,)]/[V(L)  (T,  -  Tj)]  (15) 

=  [V(IOO)  -  V(— 50)]  / [V (25)  •  150] 

=  [V(T  1)  +  V(T,)  -  2  •  V(T2)]/[V(T2)  •  2  (T3  -  T2)(T2  -  T,)]  (16) 

=  [V(100  +  V(-50)  -  2  •  V(25)]/[V(25)  •  2  (75)(75)] 

In  the  more  recent  works,  six  temperature  points  were  calculated  for  each  orientation  to  provide  data  for  the 
linear  regression  analysis.  The  six  temperature  points  were  100°C,  85°C,  75°C,  -25°C,  -10°C,  and  -50°C.  The 
results  were  used  to  calculate  TCF"’,  TCF'21  and  TCF131. 

(2)  Perturbation  Technique  of  Calculating  Temperature  Dependence  of  the  Rayleigh  Velocity 

Pertubation  theory  may  be  applied  to  the  problem  of  calculating  the  first,  second  and 
third  order  dependence  of  the  Rayleigh  velocity  V,  of  a  piezoelectric  substrate,  once  the  solution  to  the 
Rayleigh  wave  propagation  at  a  temperature  reference  T0  is  known. 

Pertubation  techniques  allow  calculations  of  small  changes  in  the  solutions  to  a  problem 
caused  by  small  changes  in  the  physical  parameters  of  the  problem,  once  the  solution  to  the  unperturbed 
problem  is  known  One  can  apply  perturbation  techniques  either  to  boundary  perturbation  such  as  mass 
loading  on  the  surface  or  to  volume  perturbations  such  as  adding  a  thin  conducting  layer  to  the  interior  of 
piezoelectric  substrates.  As  the  boundary  conditions  are  unaffected  by  changes  in  temperature  while 
material  constants  such  as  c„  are  temperature  dependent,  one  can  apply  the  volume  perturbation  theory  to 
the  problem. 

The  general  approach  to  the  problem  of  determining  the  temperature  dependence  of  V, 
will  be  as  follows.  First,  the  Rayleigh  wave  propagation  problem  will  be  solved  in  the  standard  way  in  its 
entirety  at  room  temperature,  T0.  Given  the  solution  of  the  problem  at  T0  and  the  dependence  of  the  physical 
constants  (such  as  c„)  on  temperature  at  T0,  one  will  apply  the  volume  perturbation  formula,  calculating  the 
temperature  dependence  of  V,.  The  dependence  of  V,  on  T  is  then  used  to  calculate  the  frequency 
characteristics  of  the  actual  device  given  the  thermal  expansion  coefficients  as  a  function  of  temperature.  At 
this  point,  the  frequency  temperature  dependence  of  the  substrate  as  a  function  of  the  crystal  cut  and 
direction  can  be  thoroughly  explored.  See  Appendix  A  for  a  complete  discussion  of  this  method. 


(3)  Approach  of  Sinha  and  Tiersten 


The  primary  difficulty  with  the  perturbation  technique  is  that  it  does  not  take  into 
account  the  change  of  coordinate  systems  induced  by  the  thermal  expansion  in  the  material.  This  comes 
about  because  the  material  distorts  as  temperature  changes.  Thus,  the  set  of  axis  is  to  which  the 
fundamental  elastic  constants  refer,  which  is  fixed  to  the  crystal,  is  no  longer  equivalent  to  the  axis  used  to 
calculate  V,.  This  problem  is  elegantly  solved  by  Sinha  and  Tiersten1’2.  The  first  simplification  which  occurs 
is  that  the  density  of  the  material  remains  constant  with  temperature.  Furthermore,  the  <*■ simply  become 
«i"’  =  a,"’,  as  the  length  in  this  coordinate  system  does  not  change.  The  only  difficulty  is  that  the  elastic 
constants  previously  used  no  longer  refer  to  the  proper  coordinate  system  and  must  be  rederived  from  the 
original  experimental  data.  This  procedure  has  already  been  carried  out  for  the  first  order  temperature 
derivation  of  quartz  but  has  not  deen  done  for  the  second  and  third  order  coefficients.  Nonetheless,  the 
procedure  yields  more  accurate  results  for  the  first  order  dependence.  See  Appendix  B  for  a  mathematical 
description  of  the  salient  features  of  this  technique. 

(4)  Differentiation  Method 

A  method  for  determining  the  theoretical  temperature  dependence  of  Rayleigh  Surface 
Waves  consists  of  formally  differentiating  the  wave  equation  and  boundary  conditions.  The  boundary 
conditions  and  wave  equation  must  be  true  at  all  temperatures,  placing  constraints  on  how  the  parameters 
of  the  wave  equation  may  vary.  In  this  technique,  the  derivatives  of  these  equations  with  respect  to 
temperature  are  set  to  zero  and  solved  for  the  velocity-temperature  dependence.  This  method  follows  the 
methods  used  by  Bechmann,  Ballato,  and  Lukaszek  to  compute  the  temperature  dependence  of  the 
fundamental  elastic  constants  from  frequency  data,  except  that  the  simplifying  assumptions  of  assuming 
bulk  wave  solutions  cannot  be  made.  This  method  was  later  used  by  Hauden3  to  search  for  temperature 
stable  cuts  of  quartz.  See  Appendix  C  for  a  discussion  of  this  method. 

(5)  Summary  of  Approaches 

Of  all  the  techniques  presented,  the  finite  difference  technique  satisfies  all  of  the  basic 
requirements  for  calculating  temperature  coefficients  of  delay.  Arbitrary  crystal  structures  may  be 


'“On  the  Temperature  Dependence  of  the  Velocity  of  Surface  Waves  in  Quartz,"  Sinha  and  Tiersten,  1978 
Ultrasonics  Symposium  Proceedings,  pp.  662-665. 

:,‘Temperature  Dependence  of  the  Fundamental  Elastic  Constants  of  Quartz,"  Sinha  and  Tiersten, 
Proceedings  of  the  32nd  Annual  Symposium  on  Frequency  Control.  1978,  pp.  150-153. 

'“Higher  Order  Temperature  Coefficients  of  Quartz  SAW  Oscillators,"  D.  Hauden,  M.  Michael,  J.  J.  Gagnepain, 
Frequency  Control  Symposium  (1978),  pp.  77-86. 
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investigated  once  the  density,  piezoelectric,  elastic,  and  dielectric  constants  and  their  temperature 
variation  are  known.  When  double  precision  on  the  computer  is  used  and  when  the  velocity  over  a  large 
temperature  range  is  calculated,  it  becomes  a  very  accurate  numerical  approach.  The  computational 
efficiency  of  the  finite  difference  method  is  not  as  high  as  some  of  the  other  techniques,  owing  to  the  need  for 
repeated  calculation  of  the  surface  wave  velocity. 

A  comparison  between  the  experimental  results  of  Schulz1  and  Motorola's  finite 
difference  program  adapted  from  Jones  et  al.2  as  shown  in  Figure  1 .  The  results  for  this  cut  and  many  other 
orientations  studied  have  been  found  to  be  very  good.  The  slight  shift  of  the  analytical  curve  versus 
experimental  curve  can  be  explained  by  a  crystal  misorientation. 


Figure  1.  Finite  Difference  Method  Versus  Experiment  ST-Cut  (Quartz) 


'“Surface  Acoustic  Wave  Delay  Lines  with  Small  Temperature  Coefficient,”  Schulz,  Manfred  B.,  Proc.  IEEE, 
Sept.  1970,  pp.  1361-1362. 

2"Numerical  Computation  of  Acoustic  Surface  Waves  in  Layered  Piezoelectric  Media  -  Computer  Program 
Descriptions”,  Jones,  Smith,  and  Perry,  Hughes  Aircraft  Company,  Final  Report,  Air  Force  Cambridge 
Research  Laboratories,  Contract  No.  F19628-70-C-0027. 
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Figure  2  displays  a  phase  delay  versus  temperature  plot  for  the  orientation  (YXwIt) 
56/40/27  in  quartz  obtained  from  Hauden's  graphs2.  Using  the  differentiation  approach,  he  calculated  «K'"  = 
0,  oi  =  -0.04  x  10  *  ppm/°C2,*»  =  —22.7  x  10  ‘  ppm/°C'  for  a  cut  close  to  (YXwIt)  56/40/27.  The  Finite 
Difference  approach  was  used  to  verify  his  first  order  calculation  but  is  in  disagreement  on  the  second  and 
third  order  calculations.  This  can  be  resolved  by  experiment. 


Figure  2.  Finite  Difference  Method  Calculation  for  Hauden’s  Cut  (YXwIt)  56°/40°/27° 

Tiersten’s  Perturbation  program,  while  very  difficult  to  extend  to  higher  orders,  is 
available  for  calculating  the  first  order  TCD 

Figure  3  shows  the  results  of  Tiersten’s  calculation1  and  compares  those  results  with 
the  Finite  Difference  method  as  well  as  experiment.  The  cut  used  for  comparison  is  the  AT-Cut.  The 
improvement  in  accuracy  is  substantial.  It  offers  the  significant  advantages  of  being  quite  cost  efficient,  as 
well  as  being  more  accurate  in  the  first  order  on  off-axis  cuts,  our  primary  area  of  interest.  In  the  task  of 
making  quartz  cuts  along  the  surface  of  zero  TCD,  accurate  values  for  the  first  order  TCD  prove  more  useful 
than  the  less  exact  values  for  the  first,  second,  and  third  order  TCD’s  provided  by  other  methods,  because  it 
is  a  necessary  (but  not  sufficient)  condition  for  zero  TCD  cut  to  get  a  zero  first  order.  Since  it  is  the 
dominating  term,  it  is  important  to  first  get  an  accurate  first  order  zero  TCD  locus.  Then  one  may  search  for 
the  intersection  of  this  locus  with  the  zero  second  and  third  order  terms  with  the  finite  difference  technique. 

'"Temperature  Dependence  of  the  Fundamental  Elastic  Constants  of  Quartz,"  B.  K.  Sinha  and  H.  F.  Tiersten, 
Proceedings  of  the  32nd  Annual  Symposium  on  Frequency  Control,  1978,  pp.  150-153. 

2 ’Higher  Order  Temperature  Coefficients  of  Quartz  SAW  Oscillators,"  D.  Hauden,  M.  Michael,  J.  J.  Gagnepain. 
Frequency  Control  Symposium  (1978),  pp.  77-86. 
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Figure  3.  Tiersten's  Method  Versus  Finite  Difference  and  Experiment 


(6)  Investigative  Approach 


A  necessary  but  not  sufficient  condition  to  find  a  temperature  stable  cut  of  quartz  is  that 
a^"’  =  0.  In  practice,  a  sufficient  condition  for  finding  a  zero  temperature  cut  is  that  of"1  =  of'1'  -  of'*'  =  0. 
Thus  every  zero  temperature  cut  must  be  on  the  locus  of  angles  which  satisfy  the  condition  of1" =0.  Thus  the 
first  problem  is  to  locate  accurately  such  cuts.  Both  the  finite  difference  approach  and  Tiersten's  method  are 
being  used  now  to  calculate  the  TCF. 


The  approach  used  in  this  program  is  to  identify  the  areas  (angles)  where  TCP",  TCP"  and  TCP"  are 
relatively  low  by  using  the  Finite  Difference  Approach  and  then  define  the  exact  orientation  that  has  zero 
TCP"  in  those  areas  by  the  Sinha  and  Tiersten  approach. 

In  case  the  condition  TCF'"=  TCF'"  =  TCP" = 0  cannot  be  met,  a  compromise  approach 
would  be  to  find  an  orientation  where  the  first  and  third  order  effects  tend  to  cancel  out  in  the  temperature 
range  of  interest  and  to  find  the  minimum  second  order  effect  orientation  amount  for  those  cuts.  This 
approach  would  provide  an  effective  low  TCF  cut  of  quartz  for  SAW  application. 
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(7)  Analytical  Results  on  Zero  TCF  on  Quarto 


IRF  standard  angle  definitions  (YX  wit)  PHI/THETA/PSI  for  quarto  were  used  throughout 
the  investigation'.  Consider  the  TCFs  to  be  functions  of  these  angles,  which  define  an  angular  volume  0<  PHI 
<  30°.  -90°  <  THETA  <  90°.  0  <  PSI  *£  190°, which  spans  the  space  of  possible  cuts  and  propagation 
directions.  The  set  of  points  at  which  TCF"  (PHI,  THETA,  PSI)  =  0  forms  a  surface  in  this  angular  volume. 
Likewise,  the  set  of  points  at  which  TCF12'  (PHI,  THETA,  PSI)  =  0  also  form  surfaces  in  this  angular  volume. 


Our  object  is  to  find  a  point  where  TCP"  =  TCF'2’  =  TCF1"  * 0.  If  the  surface  of  zero  TCF" 
intersected  with  the  surface  of  zero  TCF12’,  the  result  would  be  a  line  (or  a  point  if  the  two  surfaces  are  tangent 
to  each  other)  of  angular  points  on  which  TCF"  =  TCF121  -  0.  The  intersection  of  this  line  with  the  surface  on 
which  TCF'"  =  0  would  yield  a  single  point  at  which  TCF'" = TCF21 = TCF" = 0.  Neglecting  higher  order  terms, 
we  would  have  found  a  temperature  stable  cut. 

The  calculated  values  of  TCF"  versus  propagation  angles  are  shown  in  Figure  4.  The 
zero  TCF'"  is  identified  by  the  areas  where  TCF"  changes  sign. 


Using  the  Finite  Difference  approach  with  the  available  crystal  constants,  the  calculated 
results  show  that  the  zero  TCF'"  surfaces  do  not  intersect  with  the  zero  TCF"  surfaces,  based  on  the 
interpolated  results  of  the  10°  x  10°  x  10°  resolution.  It  is  not  likely  that  a  finer  resolution  will  provide 
contrary  information  because  TCF'"  and  TCF'2'  are  relatively  slow  varying  functions  as  shown  in  Figure  4. 


"Standards  on  Piezoelectric  Crystals  1949.”  Proc.  IRE  14,  Dec.  1949,  pp.  1378-1395. 
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Figure  4K.  Calculated  Values  of  TCP"1  Versus  Propagation  Angles 
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Figure  4W.  Calculated  Values  of  TCF"’  Versus  Propagation  Angles 
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Figure  411.  Calculated  Values  of  TCP"  Versus  Propagation  Angles 
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Figure  4MM.  Calculated  Values  of  TCF"’  Versus  Propagation  Angles 
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Calculations  were  performed  on  a  10°  x  10°  x  10°  grid  over  the  angular  ranges  0  <  PHI 
30°,  0  PSI  ^  180°,  and  -90°  ^  THETA  ^  90°.  These  angular  ranges,  due  to  the  symmetry  of  quartz,  include 
all  possible  angular  orientations.  These  initial  calculations  defined  the  “angular  volumes"  of  low  TCF 
orientations.  Calculations  were  then  performed  on  a  2.5°  x  2.5°  x  2.5°  grid  near  promising  orientations.  In 
this  way.  the  entire  angular  range  was  explored  and  a  large  computer-based  data  file  built.  Maps  of  first  and 
second  order  TCF’s  are  shown  in  Figure  5.  Zero  first  order  TCF  contours  are  drawn.  In  addition,  contours  of 
the  second  order  TCF  are  drawn.  In  addition,  contours  of  the  second  order  TCF  are  drawn  and  shaded  in 
areas  where  the  second  order  TCF  is  less  than  0.01  PPM/°C‘.  This  represents  a  substantial  improvement  over 
ST  quartz,  for  which  the  second  order  for  TCF  is  approximately  0.03  PPM/°C2.  Where  data  points  were 
missing,  lines  were  connected  by  interpolation. 
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Figure  5.  TCF  Map  of  SAW  Device  on  Quartz 

Despite  the  number  and  density  of  points  at  which  the  first  and  second  order  TCFs  were 
calculated,  wherever  TCF" 1  was  found  to  be  less  than  or  equal  to  zero,  we  found  TCF121  to  be  less  than  0.  The 
reason  for  this  probably  lies  in  the  lack  of  independence  of  the  crystal  constants  themselves.  Using  the 
volume  perturbation  approach  of  Auld',  one  obtains  equation  A-8  of  Appendix  A. 


'"Acoustic  Fields  and  Waves  in  Solids,”  Auld,  B.  A„  Vol.  II,  John  Wiley  &  Sons,  1973,  N.  Y„  p.  297. 
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For  small  changes  in  p  and  c1”,  one  might  expect  that  the  TCFs  are  well  correlated  with  the  temperature 
derivatives  of  the  fundamental  constants  for  Rayleigh  waves  in  quartz.  A  list  of  some  of  the  crystal  elastic 
constants"  temperature  derivatives  is  in  Table  1. 


TABLE  1.  CRYSTAL  ELASTIC  CONSTANTS’  TEMPERATURE  DERIVATIVES 


TC1" 

TC'2’ 

TC'31 

Ci, 

-48.5x10  h 

-107x10’' 

-70x1012 

Ci; 

-3000 

-3050 

-1260 

Cm 

-550 

-1150 

-750 

Cu 

101 

-48 

-590 

Cu 

-160 

-275 

-250 

C44 

-177 

-216 

-216 

Cm. 

178 

118 

21 

P 

-34.92 

-15.9 

5.3 

arn 

13.7 

6.5 

-1.9 

On 

7.5 

2.9 

-1.5 

Figure  6  shows  a  plot  of  these  values.  The  correlation  is  quite  high,  except  for  the  case  of 
C i4  and  suggests  that  the  TCFs  should  be  correlated  also.  Plots  ofTCF"’and  TCP21  versus  angles  shown  in 
Figure  4A  through  4NN  reflect  this  correlation. 

Of  course,  if  the  two  TCF  surfaces  do  not  intersect  to  form  a  line  of  zero  TCP"  and  TCP2’, 
we  must  use  a  slightly  different  approach  to  finding  temperature  stable  cuts.  The  largest  variation  in 
constants  occurs  for  Ci:.  As  we  are  interested  in  the  temperature  range  from  -50°C  to  100°C,  the  maximum 
deviation  of  T  from  T„  =  25°C  is  |  T  -Tu|  =75°C.  Calculating  the  changes  in  the  constant  Cu  for  this  maximum 
temperature  change  gives: 


1 

Ci: 

(I) 

1  A  Ci: 

I  =  2.25  x  10  1 

1 

Cr 

1  A  C,.-’1 

|  =  1.72x10 2 

1 

c,- 

I  a  c,;" 

|  =  5.32  x  10  4 

"  Higher  Order  Temperature  Coefficients  of  the  Elastic  Stiffnesses  and  Compliances  of  Alpha-Quartz", 
Bechmann,  Ballato,  Lukaszek,  Proc.  IRE,  Aug  1962,  pp.  1812-1922. 
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Thus  the  most  significant  term  is  the  first  order  term.  The  second  order  term  is  still  quite  significant  but  one 
order  of  magnitude  below  the  first,  while  the  third  order  term  is  almost  three  orders  of  magnitude  below  the 
first.  Since  the  delay  time  r  is  dependent  on  the  crystal  constants,  the  most  significant  term  should  be  TCP 
The  least  significant  term  should  be  TCP'1.  After  numerous  calculations  of  the  TCFs  for  many  orientations, 
this  appears  to  be  the  case. 

The  investigative  approach  used  has  been  to  first  locate  the  surfaces  of  zero  TCP"  (the 
most  significant  term)  with  the  finite  difference  program.  Near  these  surfaces  of  zero  TCP",  low  values  of 
TCP'’  are  sought,  using  already  calculated  results  of  the  finite  difference  programs.  Where  low  values  of 
TCP’1  have  been  found,  the  perturbation  approach  was  used  to  more  accurately  locate  the  zero  TCP" 
surface,  this  being  the  most  significant  term  in  the  total  temperature  dependence.  TCP"’s  are  then 
calculated  to  assure  that  their  effect  on  the  total  temperature  dependence  is  small.  To  date,  this  has  always 
been  found  to  be  the  case. 

(8)  Results  of  the  Investigative  Approach 

Table  2  consists  of  a  summary  of  the  results  of  using  the  investigative  approach 
described  in  paragraph  (7)  above.  There  are  three  areas  where  low  TCF  cuts  have  been  located.  These  areas 
are  centered  near  (YX  wit)  0/27/138,  (YX  wit)  7/27/135.5,  and  (YX  wit)  15/40/40.  These  orientations  have 
zero  TCP",  calculated  by  the  Sinha  and  Tiersten  approach,  with  TCP2’  and  TCP"  calculated  with  Finite 
Difference  approach.  These  areas  were  chosen  because  of  zero  TCP",  low  TCP21,  and  a  low  TCP"  which  can 
be  mostly  cancelled  out  by  the  first  order  TCF  if  the  propagation  direction  is  slightly  rotated  away  from  the 
zero  TCP"  direction.  Then  the  TCP2'  term  will  dominate  the  performance  characteristics.  The  angular 
resolution  in  these  areas  is  1°  x  1°  x  1°.  The  cuts  potentially  have  one  half  to  one  third  the  temperature 
coefficients  of  ST-Cut  quartz. 
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TABLE  2.  PROPAGATION  CHARACTERISTICS  OF  SELECTED  ORIENTATIONS 


Angles  of  ZTCFm, 

Degrees 

(S  &  T’s  Program) 

TCF'2,/°C2  (xIO ") 

Finite  Difference 

Program 

TCF,37°CJ  (x10‘°) 
Finite  Difference 

Program 

Phi 

Theta 

Psi 

6 

26 

136.31 

-1.4 

6 

27 

135.93 

-1.3 

0.67 

6 

28 

135.59 

-1.3 

0.57 

■1 

26 

135.99 

-1.5 

■ 

27 

135.64 

-1.4 

n 

28 

135.27 

-1.3 

0.65 

8 

26 

135.74 

-1.4 

0.65 

8 

27 

135.36 

-1.4 

8 

28 

134.97 

-1.3 

1 

26 

137.78 

-1.2 

0.68 

1 

27 

137.48 

-1.2 

0.65 

1 

28 

137.17 

-1.1 

0.67 

0 

26 

138.07 

-1.2 

0.67 

0 

27 

137.78 

-1.1 

0.68 

0 

28 

137.49 

-1.1 

0.62 

-1 

26 

138.37 

-1.2 

0.60 

-1 

27 

138.09 

-1.2 

0.62 

-1 

28 

137.80 

-1.1 

0.73 

14 

39 

40.195 

-1.0 

0.64 

14 

40 

40.415 

-1.0 

0.66 

14 

41 

40.64 

-1.0 

0.75 

15 

39 

39.79 

-1.0 

0.63 

15 

40 

40 

-1.0 

0.74 

15 

41 

40.23 

-1.0 

0.73 

16 

39 

39.4 

-1.0 

0.68 

16 

40 

39.605 

-1.0 

0.66 

16 


41 


39.825 


-1.1 


0.60 
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d.  Piezoelectric  Coupling  Factor 

The  piezoelectric  coupling  factor,  denoted  by  k2,  is  a  measure  of  the  coupling  efficiency  for  an 
interdigitated  electrode. 

We  can  express  k2  in  terms  of  the  interdigital  transducer's  input  conductance1  as 
G,„  (<on)  =  — jt2  a>„  M2Wt,,0°12  G,P„-i2  (COS  *7  7r)/K2  (COS  17  jr/2)  (18) 

where  =  2n  (2n  -  1)/xs.„  M  is  the  number  of  electrode  pairs,  w  is  the  length  of  the  electrodes,  x  is  the 
spatial  periodicity  of  the  array,  the  width  of  each  electrode  is  j?x/2,  t,1”'  is  the  value  of  the  effective 
permittivity  at  zero  velocity,  K  is  the  complete  elliptic  integral  of  the  first  kind,  and  Pm  are  the  Legendre 
polynomials.  G.  is  the  residue  of  1/(  |s  |«s). 

k;  is  defined  by 

k;  =  -2t,'°°'G,  (19) 

Thus  (18)  becomes 

6.W  =  (1/2)  7T2  <U„  M2  Wt,'001  k,2  P„-|2  (cos  T)  tt)/K2  (cos  17  jr/2)  (20) 

and  we  see  that  k;  is  related  to  the  input  conductance  and  hence  to  the  coupling  efficiency.  A  good 
approximation  to  k\  is  given  by 

k,2  ^  2(V„-V,)/V,  =  2AV/V  (21) 

where  V,  is  the  shorted  Rayleigh  wave  velocity  and  V0  is  the  free  Rayleigh  wave  velocity. 

It  is  seen  that  the  quantity  aWV  is  an  important  parameter  as  it  is  a  direct  measure  of  the 
coupling  efficiency  and  conductance  of  an  interdigitated  electrode  pair  on  a  piezoelectric  substrate. 
Furthermore,  aV/V  is  easily  determined  once  the  shorted  and  free  surface  problems  have  been  solved. 
Programs  are  available  for  calculating  the  Rayleigh  wave  velocity  with  either  shorted  or  free  boundary 
conditions. 

The  shorted  velocity  calculations  assumes  a  massless,  perfectly  conducting  layer  on  the 
surface  of  the  crystal.  The  boundary  condition  which  must  be  satisfied  is  that  <t> =0  at  x? =0.  The  free  surface 
boundary  conditions  are  such  that  the  potential  <*>  and  the  normal  component  of  the  displacement,  0,,  are 
continuous  at  the  surface.  Furthermore,  must  satisfy  Laplace's  equation  above  the  surface,  resulting  in 

<P  =  <t><,  exp(-kxj)  exp  (— ik(xi  -  Vt)),  (x3^0)  (22) 


"'Surface  Wave  Filters,”  Matthews,  Herbert,  John  Wiley  &  Sons,  New  York  (1977). 


In  both  calculations  the  mechanical  boundary  conditions  are  the  same,  that  there  be  no  force 
component  in  the  x.  direction,  or 

(23) 

Tii  —  Tj;  —  Tji  —  0  at  Xi  —  0 

AV/V  for  different  crystal  types  with  various  orientations  were  calculated.  Results  have  been 
found  to  be  in  good  agreement  with  experiments. 

For  example,  on  ST-quartz, 

V„  =  3.1586  x  10’ m/s 

V.  =  3.1569  x  10' m/s  (24) 

aV/V  =  5.4  x  10  J 

To  insure  the  suitability  of  the  cuts  described  in  paragraph  (8)  above  for  SAW  applications,  the 
coupling  coefficients  have  been  calculated  for  these  cuts  and  are  summarized  in  Table  3. 

TABLE  3.  PROPAGATION  CHARACTERISTICS  OF  SELECTED  ORIENTATIONS 


TABLE  3.  PROPAGATION  CHARACTERISTICS  OF  SELECTED  ORIENTATIONS  (CONT) 


Angles  of  ZTCF"',  degrees 
(S  &  T’s  program) 

Velocity 

(msec) 

K; 

(x10  ') 

. 

Power  Flow 

Angle  (Degrees) 

Phi 

Theta 

Psi 

14 

fl 

3298.60 

0.96 

-7.7 

14 

1 

3306.67 

0.96 

-8.1 

14 

41 

3315.19 

0.94 

-8.6 

15 

39 

39.79 

3301.82 

0.96 

-7.8 

15 

40 

3310.14 

0.94 

-8.3 

15 

41 

40.23 

3319.09 

0.98 

-8.6 

16 

39 

39.4 

3305.38 

0.96 

-8.0 

16 

40 

39.605 

3314.03 

0.98 

-8.4 

16 

41 

39.825 

3323.15 

0.92 

-9.0 

e.  Power  Flow  Angle 

The  power  flow  angle  for  a  particular  direction  of  propagation  is  an  important  design 
parameter.  While  the  phase  fronts  always  remain  parallel  to  the  source  transducer,  the  wave,  as  a  whole, 
does  not  propagate  perpendicular  to  the  wave  fronts  (see  Figure  7).  This  is  a  characteristic  of  anisotropic 
substrates  where  the  phase  velocity  is  asymmetric  about  the  propagation  direction;  i.e.,  V(i£  +  a<a)  *  V(<£  - 
M).  The  major  problem  which  arises  is  that  the  acoustic  beam  may  steer  off  of  the  desired  propagation 
track,  missing  the  output  trandsucer  unless  it  is  properly  designed. 


The  power  per  unit  width  carried  in  a  surface  wave  is  found  by  integrating  the  mechanical  and 
electrical  Poynting  vectors  to  obtain 


0 


s 

— oo 


T,j  mj*  dx3 


oo 


— oo 


<f>0,  dx3 


.1  =  1,2 


(25) 


P,  gives  the  power  flow  perpendicular  to  the  wave  front  and  P2  gives  the  power  flow  parallel  to 
the  wave  front.  P.  =  0  for  the  Rayleigh  wave  which  is  confined  to  the  surface.  The  power  flow  angle  may  be 
defined  as 

0  =  arctan(P:/Pi)  s  P2/Pi  for  P2  P,  (26) 


The  power  flow  angles  are  calculated  using  either  the  perturbation  programs  or  the  finite 
difference  routine.  Table  4  gives  results  of  calculations  for  the  ST-Cut  quartz.  Note  that  for  <//  =  0,  the  power 


l 
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flow  angle  is  zero,  within  the  single  precision  accuracy  used,  as  a  result  of  the  crystal  symmetry.  Power  tlow 
angles  as  high  as  20  degrees  are  not  uncommon  on  quartz. 


TABLE  4.  POWER  FLOW  ANGLES  ST-CUT 


Transverse/Incident  Power 

Power  Flow  Angle  0 

0° 

6.159  x  10  ' 

0 

10° 

6.145  x  10  ' 

3.5  degrees 

20° 

9.502  x  10 

5.4  degrees 

The  beam  steering  angle  will  be  calculated  for  the  selected  cut  of  quartz  with  TCF  =  0.  The 


analytical  results  will  be  compared  with  experimental  measurements.  Table  3  contains  the  results  of 
calculations  for  the  selected  orientations  described  in  paragraph  (8)  above. 


f.  Proximity  and  Excitation  Strength  of  Bulk  Acoustic  Waves  (BAW)  Spectrum 

The  purpose  of  this  task  is  to  analytically  determine  the  strength  of  spurious  signals  caused  by 
BAW.  With  this  knowledge,  one  can  predict  the  quality  of  the  SAW  device  which  has  been  selected  for  its  zero 
TCD  cut. 


Interdigital  transducers  do  not  couple  all  of  their  field  energy  into  surface  waves.  Bulk  waves 
are  also  generated  at  various  frequencies  in  the  piezoelectric  medium.  These  bulk  waves  can  bounce  off  of 
the  bottom  surface,  undergoing  mode  conversions  in  the  process,  and  be  received  by  the  output  transducer, 
resulting  in  unacceptable  spurious  signals.  In  most  cases,  this  problem  can  be  virtually  eliminated  by  proper 
preparation  of  the  bottom  surface. 


A  transducer  with  periodicity  a  excites  surface  waves  at  center  frequencies 


where  Vv\»  is  the  surface  wave  velocity  and  n  is  an  integer,  n =0, 1, 2 . The  coupling  to  higher  harmonics 

depends  on  transducer  design.  Bulk  waves  are  radiated  into  the  medium  at  an  angle  e  at  a  center 
fundamental  frequency 

ft.  =  V„(0)/[ x  cos(0)]  (28) 

where  Vt-(fl)  is  the  velocity  of  the  bulk  wave  in  that  particular  direction.  Equation  28  is  a  statement  of  the 
Bragg  condition  (see  Figure  8).  If  the  medium  is  isotropic  for  each  mode,  Vh  is  a  constant.  In  the  general 
case,  Vh  depends  on  the  direction  of  propagation  and  hence  on  0.  To  calculate  V„(0),  the  matrix  of  elastic 
constants  are  rotated  through  the  angle  0  in  the  standard  way.  The  quasi-longitudinal  velocity  and  quasi¬ 
shear  velocities  are  calculated  as1: 

Vi  =  i[Cu(0)  t  6n*(0)/ti i(0)]/p)'  ' 

Vs  „  =  |[C„„(0)  +  ei63(0)/t.i(0)]/p),J  (29) 

Vs\  —  |[C<<(0)  +  6i<‘(0)/tn(0)]/p}1 ' 

- - A - - 


Figure  8.  Bragg  Condition  for  Excitation  of  BAW 

The  inverse  values  of  these  velocities,  plotted  as  a  function  of  0,  form  the  inverse  velocity  curves. 

Inverse  velocity  plots  as  a  function  of  propagation  are  particularly  useful,  for  as  long  as 
Vh/cos  0,  the  effective  velocity  of  the  bulk  wave  on  the  surface,  is  greater  than  Vsa»,  the  effects  of  bulk 
mode  generation  may  be  suppressed  by  appropriate  design  and  by  proper  conditioning  of  the  bottom 
surface.  Should  Vh/ cos  0  be  less  than  VSAw  for  any  bulk  mode,  the  SAW  may  couple  to  the  BAW,  resulting  in  a 


'"Physical  Acoustics,"  Mason,  Academic  Press,  NY,  Vol.  1,  Part  A,  1964 


leaky  surface  wave.  Inverse  velocity  plots  have  been  made  for  orientations  with  promising  SAW 
temperature  characteristics  to  check  for  the  possibility  of  leaky  surface  waves  and  minimum  values  of  fh 
calculated  and  compared  with  1  SA  VS . 

The  polar  plots  of  the  inverse  velocities  for  a  (YX  wit)  0/27/137.8  and  7/27/135.59  are  shown  in 
Figure  9  and  Figure  10,  respectively.  The  values  of  (1  /vhulk)  for  the  two  shear  modes  and  one  longitudinal 
mode  for  different  propagation  directions  (0)  into  the  crystal  are  also  shown  in  Figures  9  and  10.  The 
inverse  surface  wave  velocity  for  0/27/137.8  is  3.06x  10  4,  that  for  7/27/137.8  is3.03  x  10 4.  These  values  are 
larger  than  the  maximum  (1  /  vhulk)  cos  0  (<2.9  x  10'4),  therefore,  the  analysis  indicates  that  a  leaky  mode  does 
not  exist  (see  Table  3  for  the  SAW  velocities). 


BULK  WAVE  INVERSE  VELOCITIES 

x10-3  FOR  (YXWLT)  00/  27  00/137  80 


-03  -02  -0.1  00  0.1  02  0.3 

X10‘3  (1  /SEC  r1 


Figure  9.  Polar  Plots  of  Inverse  Velocities  for  a  (YXwIt)  0/27/137.8 
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where 


ru  —  On  b2  +  2Cd  b  +  Cn 

r;;  =  C44  b2  +  2c4»  b  +  cm 

Tjj  —  Cn  b*  +  2C,s  b  +  Css 

r  12  =  C4s  b2  +  (c  14  +  cs»)  b  +  Cit 

Tn  =  C»  b2  +  (Ci)  +  Css)  b  +  Cis 
Pa  =  C)4  b2  +  (Cm  j"  C4s)  b+  Cs* 

r«4  =  -(*n  b2  +  2«,sb  +  e„)  (31) 

Tu  =  6)s  b2  +  (eis  +  6)i)  b  +  e.i 
r  24 =  eM  b2  +  (Cu  +  6m)  b  +  c  i« 

r )4 =  e))  b2  +  (fin  +  6)s)  b  +  en 
The  boundary  conditions  become,  in  matrix  form, 

...(C)).l+C)),)b'm’ai'^’,+(el))+e)))b(m')a4(",■..  !  Ci  (32) 

...(C)iii +Csii,  b«-’) «)'"’+ (e,„ + ess)  b"”)  a4,m' ...  c2 

...(C)2n+  c2)i)  b'*’)  a/*' + (e,» + emb'"*)  oV*' ...  Ci  =0 

...(e).,  +  6)1,  b1"”)  afm)-  (tji  +  e))  b,m) -  ie„)  a/1"'..^  [  C« 

The  definition  of  the  constants  and  derivation  of  the  equations  are  given  in  Appendix  C.  The 
condition  for  decoupling  requires  that  the  matrix  in  equation  30  has  zero  elements  such  that  independent, 
non-Raleigh  wave  solutions  may  exist.  Equation  30  may  decouple  in  many  ways.  If,  for  example,  r,2  =  rM  = 
r,4  =  0,  u,  and  u,  are  found  to  be  coupled  and  u2  and  <*>  are  found  to  be  coupled,  however  u,  and  u3  are 
decoupled  from  u2  and  <*>.  If  the  physical  constants  are  such  that  these  two  solutions  are  not  coupled  through 
the  boundary  conditions  (equation  32)  then  we  find  the  Rayleigh  like  wave  (u<  and  u,)  is  not  coupled  to  <t>,  the 
potential  term,  and  cannot  be  excited  by  electrodes  in  this  cut. 

The  displacement  u2  which  is  coupled  to  <t>  is  called  the  Bleustein  —  Gulyaev  wave  and  is 
excited  by  interdigital  transducers.  These  two  waves  (u3,  u,  and  u2,  <t>)  are  degenerate  as  they  propagate  with 
the  same  velocity.  The  simplest  method  to  determine  whether  equation  30  has  decoupled  for  a  particular 
orientation  is  to  calculate  the  matrix  of  equation  30.  This  is  presently  accomplished  using  the  Rayleigh  wave 
velocity  calculation  program,  which  calculates  and  prints  the  matrix  in  equation  30.  Because  of  the  variety 
of  special  cases'  which  may  arise,  some  of  which  may  ultimately  prove  useful,  each  case  in  which  we  find 
the  wave  equation  decouples  will  be  considered  on  an  individual  basis. 

'“Propagation  of  Piezoelectric  Surface  Waves  on  Cubic  and  Hexagonal  Crystals,"  Cambell  &  Jones,  J.A.P., 
Vol.  41:2796-2801  (1970). 
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The  analytical  results  of  this  section  allow  us  to  determine  which  modes  can  be  generated  on 
the  selected  zero  TCD  cut  of  quartz.  At  the  three  areas  tabulated  in  Table  2,  only  Rayleigh  wave  modes  were 
found. 

h.  Sensitivities  Due  to  Crystal  Misorientation 

In  cutting  quartz  and  aligning  masks  on  it,  there  is  always  some  maximum  achievable 
accuracy.  Thus  it  is  useful  to  know  how  all  of  the  acoustic  quantities  considered  vary  with  angle.  Quantities 
such  as  TCD,  phase  velocity,  power  flow  angle,  AV/V,  bulk  wave  spectrums,  and  bulk  wave  velocity 
surfaces,  are  of  interest  to  this  program.  These  quantities  can  be  accurately  determined  by  directly 
calculating  the  quantities  at  y>  =  (y>„  +  Ay>),  0  =  (ft,  +  se),  and  i h = (yr„+ Ayr)  with  the  same  computer  program 
discussed  in  paragraph  1,  where  y>.„  ft,  and  yr„  are  the  desired  angles,  a y>,  A0  and  a  yr  are  the  actual 
directional  deviation  from  the  desired  direction  due  to  the  fabrication  tolerance.  It  is  important  to  calculate 
the  sensitivity  of  the  parameters  to  the  crystal  misorientation;  i.e.,  amount  of  change  of  a  function  as  a  result 
of  small  angular  misorientation. 

For  case  of  the  power  flow  angle  (PFA) 

d[PFA  (4>  =  0.  9  =  42.75,  yr  =  10°)]/d*  (33) 

a  (PFA  (yr  =  20°)  -  PFA  (yr  =  0°))/20° 

a  0.27°/degree. 

The  quantities  d[PFA  (y> = 0°,  e = 42.75,  = 1 0°)] / dy>  and  d[PFA(y> = 0°,  6 = 42.75°,  >p = 1 0°)]  /do  are  computed  in 
an  equivalent  manner. 

In  the  case  of  TCD: 

dfrco)  =  +  *  ^eeid*  1341 

=  [TCD(y>„,  ft,,  yr0)  -  TCD(y>0,  ft, + Aft  yr0)]/ Ad 

de  (35) 

=  rrC0(<fro,  e0,  yr0)  -  TCD(y>0 + Ay>,  ft,,  yr„)]/ A<t> 

=  ^('^(<^>o’  _  TCD(y>0,  d0,  yr0 + Ayr)]/ Ayr 

All  the  TCD  (y>,  ft  yr)  and  TCD  (y>,  ft  yr  +  Ayr)  etc.,  are  calculated  with  the  computer  programs  discussed  in 
paragraph  1.  Therefore,  all  the  a(TCD)/3ft  a(TCD)/9yr,  a(TCD)/ay>  can  be  accurately  calculated.  Hence  the 
sensitivities  due  to  crystal  misorientation  are  determined. 

Calculation  of  the  angular  dependence  on  the  first,  second,  and  third  order  TCDs  is,  of  course, 
our  primary  task.  01  these  three  quantities,  the  first  order  TCF  is  most  sensitive  to  angular  variation.  The 
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angular  dependence  on  these  parameters  will  be  calculated  in  the  same  way  as  all  of  the  other  quantities  but 
on  a  much  smaller  angular  grid  (A0,  ±<t>  and  a< li)  about  the  zero  TCD  locus.  The  size  of  this  grid  will  vary, 
depending  on  the  magnitude  and  smoothness  of  the  variation  about  each  point  on  the  locus  of  zero  first 
order  TCD  orientations  as  verified  by  experiment  and  calculation.  In  case  the  function  is  smooth  or  linearly 
varying  with  angles,  large  A  angles  may  be  used. 

By  performing  the  above  mentioned  calculations,  we  have  precise  information  on  the 
sensitivities  due  to  crystal  misorientation.  This  information  will  allow  us  to  impose  a  practical  tolerance 
limit  on  fabrication  and  still  be  able  to  achieve  the  required  superior  performance  specification. 

Quantities  such  as  velocity  (Table  3),  power  flow  angles  (Table  3),  BAW  spectrum  (Figures  50 
and  51),  coupling  coefficients  (Table  3),  and  second  and  third  order  TCFs  (Table  2)  do  notvary  quickly  with 
angle.  This  is  not  the  case  for  TCP".  Table  5  contains  a  summary  of  aTCF^/a^.  The  large  values  of 
aTCF"7dt/r  impose  strict  fabrication  tolerances  on  the  SAW  cuts  and  mask  alignment.  Fabrication  accuracy 
to  within  6  minutes  is  required  to  keep  the  total  temperature  variation  due  to  TCF1"  within  45  ppm  for 
aTCF'"/ay»  =  3(PPM/°C)/degree  over  the  temperature  range  -50°C  to  100°C.  Table  6  contains  summaries  of 
aTCP"/9<*>  and  aTCP'790.  These  values  impose  fabrication  tolerances  on  the  rotated  quartz  plate  angles  <t> 
and  e  of  12  minutes  to  keep  the  total  temperature  variation  due  to  aTCP"(15/40/40)/a4>  within  45  ppm  over 
the  temperature  range  -50°C  to  100°C.  This  linear  temperature  variation  may  be  compensated  for  by  varying 
<b  on  any  particular  cut  if  all  other  cut  parameters  vary  slowly  with  angle. 


TABLE  5.  aTCF"’/a*  FOR  SELECTED  CUTS 


Angles  of  ZTCF"1,  degrees 
(S  &  T’s  program) 

aTCF<"/a* 

Phi 

Theta 

Psi 

6 

26 

136.31 

+2.7  (ppm /C°)/ degree 

6 

27 

135.93 

+2.7 

6 

28 

135.59 

+2.7 

7 

26 

135.99 

+2.7 

7 

27 

135.64 

+2.7 

7 

28 

135.27 

+2.7 

8 

26 

135.74 

+2.7 

8 

27 

135.36 

+2.7 

8 

28 

134.97 

+2.7 

1 

26 

137.78 

+2.8 

] 


TABLE  5.  aTCF;"/a*  FOR  SELEC1 

rED  CUTS  (CONT) 

Angles  of  ZTCF'",  degrees 
(S  &  T's  program) 

aTCF"7d* 

Phi 

Theta 

Psi 

1 

27 

137.48 

+2.8 

1 

28 

137.17 

+2.8 

0 

26 

138.07 

+3.0 

0 

27 

137.78 

+3.0 

0 

28 

137.49 

+3.0 

-1 

26 

138.37 

+3.0 

-1 

27 

138.09 

+3.0 

-1 

28 

137.80 

+3.0 

14 

39 

40.195 

-3.5 

14 

40 

40.415 

-3.5 

14 

41 

40.64 

-3.5 

15 

39 

39.79 

-3.5 

15 

40 

40 

-3.5 

15 

41 

40.23 

-3.5 

16 

39 

39.4  | 

-3.7 

16 

40 

39.605 

| 

-3.7 

16 

41 

39.825 

-3.7 

TABLE  6.  aTCF"7d<*>  AND  aTCF'Vdtf  FOR  SELECTED  CUTS 


Angles  of  ZTCF"1 
(S&T’s  Program),  Degrees 

aTCF’/a* 

aTCF"7ae 

Phi 

Theta 

Psi 

7 

■■ 

135.64 

-0.7(ppm/C°)/degree 

-0.5(ppm/C°) /degree 

0 

MM 

137.78 

-0.8 

—0.8 

15 

40 

40.00 

+1.5 

-0.7 

WAFER  FABRICATION 


During  the  developmental  phase  of  the  doubly  rotated  cut  SAW  Devices  program,  a  large  number  of 
wafers  with  30  different  crystal  orientations  will  be  fabricated  and  evaluated.  To  perform  this  task  with  good 
quality  control,  minimum  cost  and  in  a  short  time,  Motorola  has  developed  the  techniques  and  processes  to 
fabricate  the  doubly  rotated  cut  of  quartz  wafers  internally.  The  crystal  boules  are  supplied  by  Motorola. 
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Carlisle.  Pennsylvania;  the  X-ray  wafer  cutting  is  performed  at  Motorola's  Semiconductor  Croup;  polishing, 
fabrication  and  testing  are  performed  at  the  Motorola  Government  Electronics  Division.  The  methods  used  to 
fabricate  the  quartz  wafers  are  described  in  the  following  sections. 

A  computer  program  was  developed  to  calculate  the  incident  angles  and  reflected  angles  for  any 
given  cut  of  crystal.  The  basic  mathematical  relations  were  derived  in  R.A.  Heising's  "Quartz  Crystals  for 
Electrical  Circuits.”  The  relations  were  formulated  to  computer  programs.  The  flow  diagram  is  shown  in 
Figure  1 1.  The  results  are  shown  in  the  printout  (see  Appendix  F).  The  incident  angle  G.  exit  angle  (G  +  G*),  are 
defined  for  each  plate  position.  The  plate  positions  are  determined  by  the  direction  of  the  rotated  axis.  (±1, 
±2,  ±3),  relative  to  the  measuring  stage.  Once  the  angles  of  (Y X  wit)  ^/e/iA  are  defined,  all  of  the  reflection 
angles  will  be  calculated  for  the  different  reflection  planes,  defined  by  Miller  indicies.  The  useful  reflection 
planes  that  provide  intense  reflection  with  low  skewed  angles  (<5  degrees)  are  selected  for  printout.  These 
results,  after  modification  due  to  change  in  standards,  are  used  to  check  the  accuracy  of  the  cuts. 

a.  Face  Definitions 

The  following  is  an  outline  of  the  procedure  to  be  used  in  this  program  to  make  a  doubly 
rotated  cut  (YX  wit)  4>lol 0.  The  quartz  bars  have  four  lumbered  faces  with  the  minus  X  axis  marked.  The 
opposite  face  is  marked  by  coloring  it  with  a  magic  marker  (blue).  One  otherface  not  opposite  the  -X  face  is 
colored  also  with  a  magic  marker  (red).  The  red  face  is  now  defined  to  be  the  +Z  axis  as  in  Figure  12.  The 
direction  of  the  +Y  axis  may  be  found  by  using  the  right-hand  rule  (+X  crossed  into  +Y  gives  +Z).  The  +Y  axis 
will  lie  along  the  length  of  the  crystal.  Note  that  there  are  two  ways  to  set  up  the  axes  on  the  crystal 
corresponding  to  choice  of  the  red  face  (see  Figure  12). 

b.  Running  the  X-Ray  Program 

The  X-ray  orientation  program  XRAY  is  run  with  the  angles  <t>lel0.  The  sequence  of 
instructions  on  the  Honeywell  560  under  the  CP-V  operating  system  is  as  follows: 

!SET  F:103  DC/MILDAT  .538;  IN 

(This  instruction  assigns  to  unit  103  the  file  containing  all  of  the  Miller  indicies  to  be  searched.) 

!SET  F104  DC/ERRDAT  538;IN 

(This  instruction  assigns  to  unit  104  the  file  containing  all  of  the  angle  perturbations  defined  in 
Heising  for  use  with  the  Laue  photographs.) 

IXRAYL  M.538 

(This  instruction  loads  and  begins  execution  of  the  program.) 

The  program  will  ask  for  the  angles  <t>lel 0.  These  angles  are  entered  in  3G  format  [i.e.,  10.0, 3.0, 
0.0  (carriage  return)]  and  the  results  of  the  X-ray  analysis  printed. 
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Figure  11.  Flow  Diagram  of  Calculation  for  Reflected  Angles  in  Doubly  Rotated  Cut  of  Quartz 
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BLUE  FACE  ON  FRONT  SURFACE 
RED  FACE  ON  TOP  SURFACE 


.  "ilO-S 


Figure  12.  Two  Possible  Ways  of  Defining  the  Crystal  Axes  on  a  Lumbered  Bar 
c.  Analysis  of  X-Ray  Program 

The  X-ray  program  used  follows  the  development  of  Heising’s  which  uses  a  book  written 
before  the  1949  IRE  Standard  we  use  now.  His  X,  Y,  and  Z  axes  will  be  denoted  as  X,  Y,  and  Z.  The  relation  to 
our  standard  is  as  follows: 


Heising 

1949  IRE  Standard 

+X 

-X 

+Y 

-Y 

+Z 

+Z 

Furthermore,  Heising’s  incident  beam  comes  from  the  left,  while  our  incident  beam  comes  from  the  right. 
When  these  differences  are  accounted  for,  the  “position  chart"  shown  in  Figure  13  must  be  used  instead  of 
the  charts  Heising  uses  in  order  that  our  X-ray  machine  and  the  1949  standard  may  be  used.  P,,  P2,  and  P,  in 
the  “position  chart”  are  the  +X,  +Y,  and  +Z  axes  of  the  rotated  plate.  An  important  note  of  caution:  the  <*> 
crystal  face  off  of  which  we  reflect  the  X-ray  beam  has  Y  or  -Y  as  its  normal.  Since  we  may  not  assume  that 
the  crystal  faces  are  parallel,  we  must  propagate  the  surface  wave  on  this  face.  This  is  important,  because  to 
achieve  the  same  crystal  properties,  we  must  rotate  by  <1/  about  the  +Y  axis,  resulting  in  a  change  of  the  sign 
of  the  third  rotation  (see  Figure  14). 


DIRECTION  OF  MASK  ALIGNMENT 
ANGLE  ROTATION  0 


2S10-S 


Figure  T4.  Direction  of  Mask  Alignment  Angle  Rotation  ^ 


d.  Determination  of  the  Actual  <i>  Angle 

The  third  angle  of  rotation  must  be  measured  from  a  reference.  The  reference  used  is  the  face 
opposite  the  red  face  of  the  crystal  (-Z  face).  To  determine  the  actual  orientation  of  the  finished  crystal  edge, 
we  must  know  the  orientation  of  this  plane  exactly  (see  Figure  15). 


Figure  15.  Crystal  Rotation  (YX  wit)  <*>/«/ 0 
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If  the  angles  a  and  p  in  Figure  16  are  measured  using  the  (003)  plane,  which  is  parallel  to  the  X  and  Y  axes, 
is  given  by: 


=  -tan  '[(-cos  <t>  sin  a  -  sin  <*>  sin  p)l( sin  0  sin  <#>  sin  a  -  sin  0  cos  <*>  sin  p  +  cos  0  \/l  -sim’a  -  sin:/3)] 

(36) 


Thus,  given  a  <jj  desired  which  we  wish  to  obtain,  we  must  actually  rotate  the  mask  by  t/fKluii  =  t/»de.  -Ai/»  as 
shown  in  Figure  17. 

e.  Determination  of  a  and  p 

The  most  appropriate  crystal  plane  to  use  for  determining  a  and  p  is  the  003  plane.  This  plane 
has  all  of  its  atoms  lying  in  planes  perpendicular  to  the  Z  axis,  and  for  a  perfect  Z  face  the  X-ray  deflection 
angles  will  be: 

G  =  Angle  of  incidence  =  25  degrees  19  minutes 
G  +  G’  =  Exit  angle  =  50  degrees  38  minutes 

For  all  positions  on  the  X-ray  machine  with  the  +2  or  -Z  face  being  X-rayed.  To  determine  a  and  p\or  the  -Z 
face,  we  use  the  following  relationships. 

a  =  1/2  (G  •  x  toward  im  inf  vi  —  G»x  iowarI)  incident) 


P  —  1/2  (G.V  TOWARD  INCIDENT  ~  G-V  IOWARD  INCIDENT 


(37) 

(38) 


TRUE  INTERSECTION 


To  determine  a  and  p  for  the  +Z  face  we  use 

a  —  M2  (G.\  io  isiimisr  —  G-x  to  incidint)  (39) 

P  —  1/2  (G-v  io  iscims  r  —  G*y  to  im  idknt)  (40) 

a'  and  p'  are  the  a  and  p  we  would  have  measured  for  a  -Z  face  that  is  parallel  to  the  +Z  face  measured, 
which  means  that  equation  (36)  is  still  valid  for  a  and  p'. 

The  quartz  crystal  will  now  have  the  orientation  (YX  wit)  <*> I  el  -  ±<i/( a,  p)  if  we  use  the  -Z  edge 

as  a  reference  and  (YX  wit)  <t>lel  -  Ai/i(«\  P')  if  we  use  the  +Z  edge  as  a  reference.  The  computer 

X-ray  program  for  these  angular  orientations  must  be  used  if  we  are  to  get  an  accurate  measurement.  The 
integrity  of  these  measurements  of  a  and  p  is  checked  by  comparing  1 12  (G-x  +  G-x)  and  1 12  (G  v + G*v)  with 
the  theoretical  result,  Gav<  =  25  degrees  19  minutes. 

f.  Procedure  Used  to  Zero  X-Ray  Machine: 

(1)  Use  reference  quartz  plate  marked  10.1 


(2)  Set  X-ray  exit  angle  to  26.5  degrees  (marked  AT  cut  on  machine) 

(3)  Set  plate  in  X-ray  machine  with  marking  toward  you  (X-ray  bounces  off  unmarked 

surface)  and  blue  line  up  (width  direction  vertical).  (See  Figure  18). 

(4)  Adjust  with  clutch  till  X-ray  reading  occurs  at  incident  angle  13  degrees  18  minutes 

(5)  Set  plate  in  X-ray  machine  with  marking  toward  you  (X-ray  bounces  off  unmarked 

surface)  and  blue  line  down  (width  direction  vertical).  (See  Figure  19). 

(6)  Angle  of  incidence  should  occur  at  13  degrees  18  minutes  +  A  angle. 

(7)  Adjust  clutch  until  angle  of  incidence  is  at  13  degrees  18  minutes  +  1/2  A  angle. 

(8)  Check  alignment  by  measuring  in  positions  of  steps  3  and  5.  The  average  of  the  two 
readings  should  be  13  degrees  18  minutes.  Any  discrepancy  represents  error  in  zeroing  the  machine  and  may 
be  made  less  than  less  than  1  minute. 

(9)  If  the  10.1  cut  is  not  available,  the  same  procedure  may  be  used  if,  for  the  cut  to  be  used, 
the  theoretical  angles  of  incidence  are  equal  for  the  particular  crystal  plane  used  in  both  positions. 


BLUE  LINE 


Figure  18.  Position  of  Step  3, 10.1  Alignment  Wafer 


Figure  19.  Position  of  Step  5, 10.1  Alignment  Wafer  Rotated  by  180  Degrees  About  -Y  Axis 
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g.  Cutting  Techniques  for  Doubly  Rotated  Orientations 

A  slicing  machine  built  by  Meyer  and  Burger  Company,  model  TS3,  was  used  to  slice  the  quartz 
boule.  The  6  inch  diamond  impregnated  blade  is  supplied  by  Maurice  Dessau,  New  York.  The  drift  accuracy 
of  the  2  inch  cut  is  approximately  3  mils  (~5  minutes). 

The  stage  of  the  saw  has  the  capability  to  rotate  in  two  dimensions;  it  is  designed  to  make  the 
doubly  rotated  cut.  The  adjustment  accuracies  are  graded  to  0.5  minute.  The  X-ray  beam  resolution  is  about 
2  minutes. 

The  combined  accuracy  of  the  wafers  obtained  experimentally  has  been  better  than  15 
minutes.  Iterative  adjustment  of  the  cutting  stage  can  bring  the  cutting  accuracy  to  within  10  minutes.  The 
cutting  procedures  are  the  following. 

(1)  Heat  brass  plate,  glass  plate  and  quartz  boule. 

(2)  Melt  wax  on  brass  plate  and  mount  glass  plate. 

(3)  Melt  wax  on  glass  and  mount  quartz  boule  on  glass. 

(4)  Let  cool  to  room  temperature. 

(5)  Screw  mount  brass  plate  to  slice  machine. 

(6)  Initial  cut  of  boule  along  the  XZ  plane  (or  directly  to  the  desired  angle  if  the  boule  was 

lumbered). 

(7)  Determine  true  atomic  plane  using  diffractometer  and  X-ray  computer  program. 

(8)  Adjust  stage  to  correct  for  true  atomic  planes,  and  obtain  the  desired  rotation. 

(9)  Perform  doubly  rotated  cut  of  quartz. 

(10)  Measure  the  reflection  angle  and  check  against  the  calculated  result  of  the  X-ray 
program.  Iterative  adjustment  of  stage  can  be  done  if  cut  accuracy  is  critical. 

These  procedures  are  presently  used  to  cut  the  wafers  with  low  TCF  orientation  required  for 

this  program. 

h.  Wafer  Polishing 

The  polishing  process  for  the  quartz  wafers  was  established  during  this  period.  The  procedure 
includes  the  following: 

(1)  The  crystal  axis  orientation  is  marked  on  the  back  surface  of  the  wafers  with  marking  ink 
and  baked  dry.  The  markings  have  been  demonstrated  to  stay  on  the  crystal  throughout  the 
photolithographic  process  for  easy  identification. 
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(2)  The  corners  of  the  wafers  are  ground  prior  to  lapping  and  polishing  to  prevent  corner 

breakage. 

(3)  The  polishing  is  a  two  step  process.  The  first  step  is  a  15  lapping  to  obtain  surface 
flatness  with  Microgrit,  type  WCA,  Size  15,  on  a  LAPMASTER,  Model  24.  The  second  step  is  to  polish  the 
surface  with  Cerium  oxide  on  a  LogiTech,  Ltd.  Model  PM2  polisher.  The  polish  wheel  is  made  of  pitch  wax 
supplied  by  Hacker  Instruments,  Inc.  The  polished  wafers  have  a  surface  finish  with  no  observable  grains  or 
pits  under  a  50X  microscope.  The  wafers  are  suitable  for  surface  wave  application.  X-ray  orientation  was 
performed  before  and  after  the  lapping  and  polishing  processes.  The  change  is  normally  less  than  6  minutes. 
At  the  present  time,  wafers  of  the  following  cuts  are  polished;  these  are  <t>lehii  of  6/27/135.9, 6/28/135.6, 
7/27/135.6.  7/28/135.3,  8/26/135.7,  8/27/135.4,  and  8/28/134.9. 

In  the  coming  period,  SAW  devices  will  be  fabricated  on  these  wafers.  Experimental 
results  will  be  obtained  and  compared  with  the  theoretical  calculations. 
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SECTION  III 

CONCLUSION 

The  temperature  coefficients  of  frequency  (TCF)  have  been  analyzed  for  doubly  rotated  cuts  of 
quartz  for  surface  acoustic  wave  devices.  The  analysis  procedure  is  shown  below: 

(1)  Using  Finite  Difference  program,  identify  angular  areas  with  zero  TCF1"  and  low  TCF'21  on 
10°  x  10°  x  10°  grid. 

(2)  Calculate  in  these  areas  with  2.5°  x  2.5°  x  2.5°  grid  to  locate  the  rotation  with  minimum 
TCF0’  which  has  TCF"’ =  0. 

(3)  Identify  the  accurate  zero  TCF("  with  Sinha  and  Tiersten's  program  in  areas  obtained 

in  (2). 

(4)  Calculate  TCF"1,  TCP21  and  TCF'31  in  selected  areas  obtained  in  (3)  with  1°  x  1°  x  1°  span. 

(5)  Calculate  the  coupling  coefficients,  propagation  velocities,  power  flow  angles  and 
inverse  bulk  wave  velocities  for  these  angles. 

The  conclusions  of  the  initial  study  are: 

(a)  The  zero  TCP"  surface  and  zero  TCP21  surface  data  does  not  show  an  intersection 
in  the  10°  x  10°  x  10°  grid.  However,  orientations  with  frequency  variation  ■»  100  ppm  from  -50°C  to  100°C, 
better  than  ST  Cut  by  a  factor  of  two,  have  been  found  in  several  areas  where  TCF"1 = 0  and  TCP2'  is  less  than 
1.5  x  10  8/C°  2. 

(b)  The  aTCF/a  angle  is  normally  very  large  in  the  areas  with  zero  TCP"  and  low 
TCP2',  therefore  accurate  crystal  orientation  is  critical.  It  is  estimated  that  the  orientation  accuracy  of  6 
minutes  is  required  to  obtain  a  TCF" 1  of  «  50  ppm  in  the  -50°C  to  100°C  temperature  range  in  the  three  areas 
investigated.  These  are  <t>lehhof  6/27/135.93  (4.45  xIO  2ppm/°C-  minute),  14/40/40.415 (5.83 x10"2ppm/°C 
■  minute)  and  0/27/137.78  (5.00  x  10~2  ppm/C°  •  minute).  Orientations  with  less  sensitivity  to  orientation 
accuracy  have  been  evaluated.  These  are  the  cuts  that  have  a[TCP"]/a</»  =  0  when  TCP"  =  0.  The  cuts 
evaluated  were  <f>/e/t/f  =  20/30/155,  20/20/150,  10/40/168.  The  TCF'2'  of  these  cuts  are  approximately 
4  x  10'7C°,  comparable  or  higher  than  ST  cut. 

(c)  Coupling  coefficients,  velocities  and  beam  steering  angles  were  defined  for  the 
selected  areas,  no  leaky  modes  were  found  in  these  areas. 

(d)  Methods  to  X-ray  orient  the  quartz  crystals,  cut  the  doubly  rotated  wafers  and 
polish  such  crystal  wafers  were  developed.  The  wafers  are  ready  for  SAW  device  fabrication.  The  effort  in 
the  coming  period  is  to  design  and  fabricate  single  mode  oscillators  and  to  experimentally  measure  the 
TCF’s  of  these  selected  orientations. 
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APPENDIX  A 


VOLUME  PERTURBATION  OF  AULD 

Perturbation  techniques,  as  used  here,  allow  calculations  of  small  changes  in  the  solutions  to  a 
problem  caused  by  small  changes  in  the  physical  parameters  of  the  problem,  once  the  solution  to  the 
unperturbed  problem  is  known. 

Our  general  approach  to  the  problem  of  determining  the  temperature  dependence  of  V,  will  be  as 
follows.  First,  the  Rayleigh  wave  propagation  problem  will  be  solved  in  the  standard  way  in  its  entirety  at 
room  temperature,  T<>.  Given  the  solution  of  problem  at  To  and  the  dependence  of  the  physical  constants 
(such  as  c„)  on  temperature  at  T0,  we  will  apply  the  volume  perturbation  formula,  calculating  the 
temperature  dependence  of  Vs.  The  dependence  of  V,  on  T  is  then  used  to  calculate  the  frequency 
characteristics  of  the  actual  device  given  the  thermal  expansion  coefficients  as  a  function  of  temperature.  At 
this  point,  the  frequency  temperature  dependence  of  the  substrate  as  a  function  of  crystal  cut  and  direction 
can  be  thoroughly  explored. 

The  Volume  Perturbed  Formula. 

Denoting  u,  (T»)  by  u„  etc.,  the  volume  perturbation  formula  is  given  by1: 

=  [Ajdy/J**  (-u*T  -u'T* 

O  O 

+  <t> *  (icuD')  +  4>  (iwD)*)  •  z  dy  (A-1) 

[A]  =  [AP  u*  •  u'  +  T*:  (ase:  T  +  Ad  •  E') 

+  E*  •  (a*t  •  E’  +  Ad:T')] 


This  equation  is  exact,  but  involves  knowing  the  solution  u'.  However,  we  can  set  u  =  u'  if  the  temperature 
dependence  of  the  solution  is  small  and  use  P  =  (1/2)  Re  J  (-u*T  +  E  H*)  dy  to  obtain  the  approximate 
solution 

A/3*  =  (w/4P)  [Ap  |  u  |  1  +  S:  Ac:S  ^  2) 

+  E*  •  At1  ■  E  +  E*  •  Ae:S  +  S:Ae  ■  E]  dy 
The  u’s  and  E’s,  and  p  come  directly  from  the  computer  solution  at  T0. 

'Auld,  B.A.,  "Acoustic  Fields  and  Waves  in  Solids,"  Vol.  II. 
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The  A  terms  can  be  written  as  expansions  of  the  form. 


A0„  =  0R  (T„)(a^"  dT  +  af'  dT2  +  o^'3'  dT3) 

(A-3) 

Ap  =  p(T0)(ai;  dT  +  a'p  dT2  +  a”'  dT3) 

(A-4) 

=  p"’  dT  +  p12’  dT2  +  p'3’  dT3 

(A-5) 

Ac„ki  =  Cijki"’  dT  +  c,)kl'2’  dT2  +  Cun'3'  dT3,  etc. 

(A-6) 

to  obtain 

A0r  =  /Jr  (T„)  ( dT  +  |3r  (To)  a/*  dT2  +  0R  (Td)  «j*  dT3) 

=  (cudT/4P)  X"  (p  (To)  a';’  |  u  | 2  +  S:c"':S  +  E*  • •  E 

O 

+  E*  •  e"’  •  S  +  S:e‘"  •  E)  dy 

(A-7) 

+  MT2/4P)  X’  (p(To)  a p]  1  u  |  2  +  S:c'2,:S  +  ...)  dy 

+  MT3/4P)  X’  (p(To)  ap'3’  |  u  |2  +  S:c,3,:S  +  ...)  dy 

O 

or 

V  =  (<o/4P/3r  (To))  f  (p(To)  I  u  | 2  +  S:cg,:S  +  ...)  dy 

(A-8) 

For  the  problem  of  quartz,  the  electrostatic  coupling  is  small  and  the  electric  terms  can  be  ignored. 
Thus  Ad  and  Ae  in  (A-8)  may  be  set  to  zero  resulting  in  the  simplified  equation. 

<y  =  (o>/4P/3r)  J*  (p(T0)  ««  |  u  | 1  +  S  :  cul:  S)  dy  (A-9) 

o 

The  integrals  are  calculable,  for  instance,  from  (14), 

X*  p  (T**)  “p  I  u  1 2  dy 

=  p(T.  )  a1;  5  X*  [  J  Cm  0|,ml  expH  0.  b""'  y)] 

[2  c:  exp  (i^b1"1  y)*]  dy 


(A-10) 


The  problem  can  easily  be  extended  to  multi-layered  media  by  simply  performing  the  integral  over 
each  layer  separately.  For  the  double  layer  (A-9)  becomes: 

Of  oo 

J  *(T.)<  |u|Jdv  +  J*  pi  (T0)  arpj  I  u  | J  dy  +  •••  (A-11) 

o  a 

Because  the  dependence  of  u,  is  given  explicitly  by  (13),  it  is  not  necessary  to  use  numerical 
integration  procedures,  as  demonstrated  in  (A-10). 

Once  the  have  been  calculated  from  (A-9)  the  av01  are  simply  determined.  Using  p  -  wl V.,  and 
(A-9)  we  determine  on.0’. 
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APPENDIX  B 

VOLUME  PERTURBATION  FORMULA  OF  SINHA  ANO  TIERSTEN 

The  approach  used  by  Sinha  and  Tiersten  includes  the  effects  of  distortion  caused  by  heating  a 
substrate  of  quartz  and  is  based  in  the  coordinate  system  to  which  the  fundamental  elastic  constants  refer. 
In  this  reference  system,  the  density  is  constant  and  the  ai'"  are  equivalent  to  the  av"'. 

In  this  new  reference  system,  the  perturbation  formula  becomes1 


(Av,/v,)  =  (1/2/?„V)  (H/N2),  (B-1) 

4  4 

N:  =  (p°rri//JR2)  £  2  (C'"”  a"”'  C,n'*  a,n'*)/(/Jro  -  0.*)  (B-2) 

m  =  l  n  =  I 

r°°  r  7T/0R 

H  —  “J  dy  J  dz  (K iz  Utt<  ^  Kry  lly,I  +  Kix  lilt  +  Kyi  Ui,y  +  Kyy  Uy,y  +  Ky*  Ui,y)  (B“3) 

0  — 7r/  /3  R 

Kl7  *  (CL7M0f  +  ACL-yMtt)  Uof,  M  (B~4) 

For  the  first  order  perturbation  in  T  we  have: 

^l7M«  ~  (^L7M>;AB  aAB  +  ^L7KM  aVK  +  ClKMV  °7k)  ^  (B-5) 

^Ct7MOf  =  W  cLyMO/dT)  dT  (B-6) 


Where  ci-km a  are  the  second  order  elastic  constants  previously  denoted  simply  as  cLyMa  and 

Ci  yMKAB  are  the  third  order  elastic  constants.  The  terms  acl7m«  as  calculated  by  Ttersfen2  is  available  only 

\  2 

to  first  order  in  T,  and  the  higher  order  elastic  constants  Cl7mvAbcd  have  never  been  determines 


'“On  The  Temperature  Dependence  Of  The  Velocity  Of  Surface  Waves  On  Quartz,”  B.K.  Sinha  and  H.F. 
Tiersten,  1978  Ultrasoncs  Symposium  Proceeding,  IEEE,  pp.  662-665. 

^‘Temperature  dependence  of  the  Fundamental  Elastic  Constants  of  Quartz,"  B.K.  Sinha  and  H.F.  Tiersten, 
Proceedings  of  the  32nd  Annual  Symposium  on  Frequency  Control,  1978,  pp.  150-153. 
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APPENDIX  C 


THE  DIFFERENTIATION  METHOD 


A  method  for  determining  the  theoretical  temperature  dependence  of  Rayleigh  Surface  Waves 
consists  of  formally  differentiating  the  wave  equation  and  boundary  conditions.  The  boundary  conditions 
and  wave  equation  must  be  true  at  all  temperatures,  placing  restraints  on  how  the  parameters  of  the  wave 
equation  may  vary.  In  this  technique,  the  derivatives  of  these  equations  with  respect  to  temperature  are  set 
to  zero  and  solved  for  the  velocity  temperature  dependence.  This  method  follows  the  methods  used  by 
Bechmann,  Ballato,  and  Lukaszek'  to  compute  the  temperature  dependence  of  the  fundamental  elastic 
constants  from  frequency  data,  except  that  the  simplifying  assumptions  of  assuming  bulk  wave  solutions 
cannot  be  made.  This  method  was  later  used  by  Hauden2  to  search  for  temperature  stable  cuts  of  quartz. 

Christoffel’s  wave  equation  can  be  written  in  matrix  form  as 


>..-pV.2 

r,2 

r„ 

r,4' 

ai 

r,2 

r22  —  pV,2 

r23 

r24 

<*2 

r,3 

r23 

1 

< 

m 

rM 

r.4 

r24 

rM 

Vti. 

04 

where 


r  1 1  —  Css  b2  +  2cis  b  +  Cu 

r22  =  C44  b2  +  2c46  b  +  Cm 
rj) =  Cu  b2  +  2C)S  b  +  Css 

rt2  -  C45  b2  +  (Cu  +  C56)  b  +  C16 

fu  =  Cjs  b“  +  (Cn  +  Css)  b  +  C15 


r2J  =  C34  b2  +  (Cj«  +  C4s)  b+  Cm 
r44=-(t3J  b2 +  2ci3  b  +  en) 

r,4  =  e3s  b2  +  (e.s  +  e3l)  b  +  e., 

r24  =  eM  b2  +  (e>4  +  e36)  b  +  e.« 

r34  =  e33  b2  +  (e.3  +  e3s)  b  +  e1} 


(C*1) 


(C-2) 


'“Higher  Order  Temperature  Coefficients  of  the  Elastic  Stiffnesses  and  Compliances  of  Alpha  Quartz,” 
Bechmann,  Ballato,  and  Lukaszek,  Proc.  IRE,  Aug  1962,  pp.  1812-1822. 

2"Higher  Order  Temperature  Coefficients  of  Quartz  SAW  Ocsillators,"  D.  Handen,  M.  Michael,  J.J.  Gagnepain, 
Proc.  Frequency  Control  Symposium  (1978),  pp.  77-86. 
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following  the  notation  used  in  the  previous  section  and  in  Matthews'.  This  equation  holds  for  each  of  the  four 
modes,  hereafter  designated  by  a  superscript  m. 

The  boundary  conditions  become,  in  matrix  form, 


+  c„„  btml)  a,1"”  +  (e,„  +  e333  b,n”)  a;1"’ ... 

...(Cjiii  +  Cjiij  b'"’)  ar  +  (em  +  e3)i  b,m’)  a,{n" ... 
. . . (C 32» i  +  C:„,  b,m))  a,1"”  +  (e.32  +  e3}2b'm’)  a4,n"  ... 

,..(e3„  +  e3l3  blm') <«3>  +  *33  b'"’  -  i«.)  a;m'... 


c; 

C2 

C3 

C4 


=  0 


(C-3) 


where  only  the  m’th  column  is  shown  and  the  Cm’s  are  the  amplitudes  of  each  mode.  This  formulation 
assumes  a  normalization  of  the  a’s 

(a,1"')2  +  (tt2ln”):  +  (a3'ml)2  +  (a4,m')2  =  1.  (C-4) 


The  first  condition  we  can  place  on  the  wave  equation  is  that  the  determinant  of  Christoffel’s 
equation  vanishes  for  all  temperature,  or 


r.i  -  pV,2 

r,2 

r,3 

r,4 

r,j 

r22  -  pV,2 

r23 

r24 

r,3 

r23 

r33  -  pV,2 

r,4 

r,4 

r24 

r34 

r  44 

For  each  of  the  four  modes,  this  equation  is  valid,  resulting  in  four  equations  in  five  unknowns,  dVs/dT,  and 
dbim’/dT  of  the  form  F,,ml  (dV,/dT,  dblm,/dT)  =  0,  m  =  1,  2,  3,  4. 

For  each  mode,  Christoffel’s  matrix  equation  (C-1)  must  vanish,  yieldinn  the  set  of  equations 

(d/dT)  (oi,ml  (rn'm)  -  pV,2)  +  «2""'  r,2,m'  +  a3,m)  ^l3(m,  +  a/1"1  riT )  =  0 
(d/dT)  (a,'"”  r,T+  a2,m’  (riT*-  pV,2)  +  a3(ml  r23<”'  +  a/"1  r(2T' )  =  0 

(d/dT)  (a,'m>  r,J'm'  +  a2,m'  r23'm’  +  a3'm'  (rM11”'  -  pV.2)  +  a/"’  ri:’ )  =  0  (C-6) 

(d/dT)  (a,""'  r,4,m'  +  aj'"”  r24,m’  +  a3'ra'  r34,ml  +  a4,ml  r44,n”)  =  0. 

This  results  in  16  equations  (four  per  mode)  and  an  additional  16  unknowns  (the  a,|m's)  of  the  form 
F:,lm|  (dV./dT,  dblm7dT,  d  a,lm7dT,  d  a2lm7dT,  d  a,'m7dT,  d  a4lm7dT)  =  0 


'“Surface  Wave  Filters,”  H.  Matthews,  John  Wiley  and  Sons,  New  York  (1977). 
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Using  (C-3)  we  obtain  the  single  equation 


(d/dt) 


(Cn,i  +  c33l3b,m>) (e,„  +  ejub""') 

(c„„  +  c3„3b"*’)  a,lml  +  (e.)i  +  e3„'B»)  a;0" 


=  0 


(C-7) 


of  the  form  F.'m'  (db"7dT,  db,:7dT . da,lm7dT)  =  0 

and  four  equations  of  the  form 


(d/dT)  £  Cm  [(c33„  +  c3Ji3b"”>)  o,tB>  +  (e,33  +  e33Jb,m))  a;-"'  ]  =  0 
(d/dT)  X  Cm  [(c3l„  +  c3i,3b(m’) +  (e.3.  +  e331b,n')  ]  =  0 

(C-8) 

of  the  form  F (db'"’7dT,  d«,lm7dT,  dCm/dT)  =  0 

(C-9) 

From  (5)  a,"*'  (da,""7dT)  +  a2|B”  (da2,mVdT)  +  a3"">  (do3'"7dT) 

(C-10) 

+  a/"”  =  0 

or 

F,,ml  (d  «i'm7dT)  =  0 

(C-11) 

Combining  the  above  results  gives  25  equations  in  the  25  unknowns. 
dV,/dT,  db"”7dT,  da,""7dT,  dCJdT 

which  are  solved  simultaneously.  Once  the  temperature  dependence  of  the  Rayleigh  wave  velocity  is  found, 
the  frequency  dependence  is  found  in  the  usual  way. 

T o  obtain  the  second  order  dependence  of  V„,  the  25  equations  are  differentiated  again.  The  values  of 
dVJdT,  db""7dT,  da,lm7dT  and  dCm/dT  previously  obtained  are  used  to  obtain  d2Vs/dT2,  etc. 
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APPENDIX  D 


SOLUTION  METHOD 

Calculating  the  velocity  of  a  Rayleigh  wave  requires  a  complete  solution  of  the  problem  to  be 
performed.  All  of  the  constants  in  equation  (13)  must  be  evaluated.  The  standard  used  at  Motorola  as 
developed  by  Campbell  and  Jones  is  outlined  below. 

First  the  fundamental  constants  are  rotated  into  the  coordinate  system  of  interest. 

Next  a  value  of  V,  is  picked.  The  coefficients  r„  of  equation  (C-2)  are  evaluated  and  the  determinant  of 
the  matrix  in  equation  (C-1)  is  set  to  zero,  as  it  must  be  if  a  solution  of  (C-1)  is  to  be  found.  This  results  in  an 
eighth  order  equation  in  b.  This  equation  is  solved  for  the  eight  complex  roots.  The  four  roots  blml  in  the  lower 
complex  plane  are  retained,  the  four  discarded  roots  not  satisfying  the  boundary  conditions  at  infinity. 
Equation  (C-1)  is  then  solved  tor  the  tour  eigenvectors  ar’.The  b,*"”  and  the  corresponding  eigenvector  a,""' 
are  substituted  into  equation  (C-3)  and  the  determinant  of  the  matrix  in  equation  (C-3)  is  evaluated.  This 
determinant  must  be  zero  for  a  solution  to  (C-3)  exist.  If  it  is  not  zero,  V,  is  varied,  and  the  whole  procedure 
repeated,  until  it  is.  Once  a  value  of  Vs  is  found  such  that  the  determinant  in  equation  (C-3)  vanishes,  the 
solution  to  (C-3)  is  found,  giving  the  values  for  Cm.  These  constants  completely  describe  the  solution  of 
equation  (13),  as  well  as  providing  the  exact  velocity. 
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APPENDIX  E 


EQUIVALENCE  OF  TEMPERATURE  COEFFICIENTS  OF  FREQUENCY  AND  DELAY 

In  the  following  section,  the  relations  relating  the  first,  second,  and  third  order  temperature 
coeffieients  of  delay  and  frequency  will  be  derived.  These  relations  show  implicitly  the  equivalence  of  each 
representation  of  the  device  temperature  characteristics,  and  justify  their  interchangeable  usage. 


In  the  text,  we  have  used  the  following  notation: 

t  =  delay  time  of  delay  line  oscillator  (E-1) 

ith  order  temperature  coefficient  of  delay  =  TCD1"  =  (E-2) 

F  =  frequency  of  delay  line  oscillator  (E-3) 

ith  order  temperature  coefficient  of  frequency  =  of"1  (E-4) 

r  =  r„  (1  +  a"'  (T  -  To)  +  aa'T( T  -  T„)2  +  a  '«  (T  -  T.)  3  +  ...)  (E-5) 

F  =  F.  (1  +  aFm  (T  -  To)  +  ap‘2’  (T  -  To)2  +  ap'31  (T  -  To) 3  +  ...)  (E-6) 


where  T  is  temperature  and  T«  is  a  reference  temperature,  25°C  in  our  case.  For  a  SAW  oscillator,  Fr  = 
constant  or 

r/r„  =  F„/F 

=  1/(1  +  or'"  (T  -  To)  +  aF121  (T  -  To)2  +  aF15'  (T  -  T0)3  +  ...)  (E-7) 

Using  the  relation 

1/(1 +X)s  1 -X  +  X2- X3  for  X«1  (E-8) 

we  can  write 

rho  =  1  -  [aF*’  (T  -  To)  +  ap'2'  (T  -  To)2  +  ap'3’  (T  -  T0)3] 

+  ((aF*11  (T  -  To)2  +  2«f‘"  (T  -  To)  ap'2'  (T  -  To)2] 

-  (of"'  (T  -  To)]3  +  higher  order  terms  (E-9) 

=  1  -  ap"’  (T  -  To)  +  (-ap12'  +  (ap'")2)  (T  -  T0)2  +  (-of'3’  +  2of‘"  ap'2'  -  W’)3)  (T  -  T.)3  (E-10) 
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r 

r  ■ 

Equating  powers  ot  T  -  T„  in  (E-5)  and  (E-10),  we  obtain 


=  -aFin 

(E-11) 

=W21  +  (ok'")2 

(E-12) 

=-of111  +  2aFl(l  of12'  -  (aF1")3 

(E-13) 

As  the  only  assumption  on  F  and  r  used  is  that  of  equation  (E-7)  which  is  symmetric  in  F  and  r,  we 
immediately  obtain 


_  <  1 )  _  (It 

Q\ - OtT 

(E-14) 

aF*2'  =  -a/21  +  (a/1)2 

(E-15) 

of1'1  =  -arl!l  +  2ar'"  ar121  -  (ar"')! 

(E-16) 

Using  (E-11)  through  (E-16),  we  can  always  relate  one  set  of  temperature  coefficients  to  the  other. 
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APPENDIX  F 
X-RAY  RUN 

The  following  is  an  example  of  the  computer  results  from  the  FORTRAN  X-ray  orientation  computer 
program  for  the  doubly  rotated  cut  7/27/0.  File  MILDAT  contains  only  the  two  Miller  indices  shown  on  the 
printout  on  this  run. 


009000* I Ail Cmiim!  C*»CC«*l«(U)V4  AilCmSiml 
<*ml**  *C)  Wumtlf  ilxl> 3*194  0  C0T4  9 
(*miM  *  ) 

(*NIn  9i)  0Mt*0t»19mi0l9mI  40  319N0  ■  9 
OVMMON  CtKV  C  *^«9liC3A  till  (  >  *  N0tiIC04  ml  3i«% 


} 

t 

APPENDIX  6 
X-RAY  PR06RAM 

The  following  is  the  FORTRAN  computer  program  used  to  calculate  X-ray  diffraction  angles  for  a 
doubly  rotated  cut  of  quartz. 


500  COMmo*  /!NT£N/REP 

ooo  integer  m 


\ 


i 

•j 

1 


■I 

j 

1 

i 


f 

S 


m «»ivm  «  * 

hi  «•%»«« 

t-  «•  ~mmi  « 

<  «•  •• 

►  x  •  %  »«u 

s  »&«5o5  <•*•««••( 

h  rv xx u«u • u i 

|?***SSl££M' 


KIM 
MU 
•  « 
(MM 


e 

M 


O  • 
I  *■  * 


•  eeee  ce  < 

■b  X  X  *  X  •.  *  •. 
'MKI9  if-  «  r*  •  < 


o  c  c  e  c  pe  c 
•occeeoc 

I,  «  I. 

O  Ml  «  |fl  «*»  _  _ 

~  -  —  X  MX  X  MM  X 


e  to  c 

.  %  ‘  V  O^MMtVoV'oV 

- Nkmt - “  * - 


t  Ml  Ml  Ml  •**  MM  Ml  1*1 


II.MI  CNO 


X  C  « 

•  m  HI 

X  C  4 

^  M  <»  am 


X  >  « 
iKM 

5«i 

M  ->  • 


X  «•  M 

*  V  ** 

♦  «*■*  W  « 

X  J  N  % 

w  W  £  A 

►  ►  X  •* 

4  4  «w| 


•* 

•- 

w 

4 

~  O  M 

M 

O 

O  4  •• 

X 

-4 

w  •  m 

X 

X 

m  4** 

• 

4 

0  w  • 

a 

«  N  4 

*» 

9 

4  4  X 

* 

m 

e  o  • 

V 

•  X  o 

w 

• 

444 

X 

4*  X  4 

•J« 

♦ 

fv4 

X  X  4 

f  r!?tt  ? « 

w 

<M%| 

2  x  x 
o  c  o 

ill: 

iuuz 


►»  x  x  «  ■  ■  hnwm  i  ► 

:  wtzxioAaA  x  •-  ~  m  mmm^m^mmm 

•  •fill-.  J  C  w5  •^%oSw?rSC«>of  X 

CKtl  141  Pm—  M  w  w  •  m  44  Z 

:ioco|mif!>ic  ■  i  o  x  o  x  o  zoo  •  Jk  1 1  mo 

.  « Nuggft Nftit no 4404  oa o«4uo mm ► 


4#«C44  4«C«44  4ecC4444t44  •••«4444Ce«  4 XC  44  4444 
444C44  4ea04«*4e«e4  44  04«4  44  4ee44404e4  44  440444 
•  •••.  «i  ft.  4  «.  «»  «  •.  4  •«•.%•.  «•.«.%  ft.  •.*•.%•.  4  «.  4  **.«.«.•  «.*•»•.  4  •.  4  4  • 

«*4ii4«i44^4»  c  —  —  mm  mm  mm  m  —  m  mm  m  mm  mm  © - • m  4k444«m 

V44444444X  MM  r  ir  mmmmmm 


4  4«r444tO*44«V  «4444  ^41^4  944  «f  4  *>4l 
•  4  4  Z4«<«r4f44irrtf  #444  44  4  441  4444M 


ammu 


ft  Ik 
to  -# 
►  Ik 
ft 

x  *- 

X  X 
U  J 
•* 

X  Ik 

* 

4*1 
•  X 
l»*-“ 


«  v 

to*  *2 
3«  Off 
-!•*•  Z> 
«  X 

>  tot# 


c 

o 


X  * 


X  Ml  •  •  M 
►  Wwo  •  O  w 


3M 

XZX 


X  X 

ft  ft 

X  to 


to  «J 

ft  U  ft  X 
HXQ 


3  ft  •  ►  ©ft  • 

u  •  ►*  to  ft  x  x. 

mV  %•  O  A  O 

ft  ft  ft  •  ~  •m  ►  ft 

k  O  ^  OO  ^  ft  X  to 

k  U  «  ••»  «  ft  k 

toft  to  x  kO«  e  zx jjw 
^*ft=»«ft  jm  v  in  m  9  •  to  « 

n  3  O  O  J  •  aft  AlMk^i  | 

w  U  •  to  *^«»toto  •  ft  ft  ft  ^ 

^  wuXkenf>k  s  k  *  •  *s  s 


X  ^ 

•  ft 

ft  ft 

•  ft 

X 


x  ©  to 

•  ft  > 

•  to  ft 
«  ii 
«  to 
«JO 

•  •k 
«ft 

•  to  k- 

•  ft 
to  •  k» 
to  o 

«  to 

U: 

«  -J  to 
to  to 

•  X 

•  to-  O 
to  to  at 

toX  to 

•  to 

to  to* 

•  x  m 

to  O  ft* 

•  5  B 

«  •>  v 
«  u 
too  o 

tog  _i 
«ft  u 
toft* 

«  to  ft 

«|»o 
*  •©  ►* 
i  «  k 

•  toft  to 
i  «  to  fr- 

*  J  o 


w 

X 


_  ss 

ft  ft  WW 

-*•-51 

iSZMsti 


u  * 

-Ift-  i 


ft  o 

VkC 
to  to  ft 


o  Q  ft  ft 
to  9  to  to 
*1  I  k 
Z>Oft-ft 


to  c  •  ft  «  m  +*  •-  «w  w  w  ••  in  to  >  ■  to  to -j  ®  © 
jMOZ«4ftkak4tiiftwft'^Vkk|  nj j 

kBk&Mk^vkZZft'B.altoSf  kktoMi  -J  J 

totoft-itototototoOtotoftftfttototoftMfte  c  to  to 

to  S  M  to  H  to  MM  to  toft  N  MM  MB  B  ft  Ml  O  Bft  ft  ft  ouu 

to  to 


kfi 

-I  to 

to  ft 


OH 

to  ft 

►  s 


is:- 

ft*  to  to 
kkf 

ssi 


>•«  c  ©  •  o  to  otootoooooototo< 


to  to  c  9tototo«c«ftfttoeeoo©ec  ooototo 
*.  to  to  ••  to  «.  «.*«•%••  •.  to  «.  x  to  «.  «.  •  a  x  < 
m  to  m  ckftftftto«*«totoMtotoitom  «k  ft  ft 
«  cfttotoftfttototoftftftft  k ttktceoacc 


V  ©to  <  ft  «  ~tof»toftftto*  ft  Cftiini 

ft  «  «  toftftftftftftftftftft  ft  e  •  o  ©  < 


‘•^SSS 


It  Mt'll 


31/lliii  ttt'lil  •  Ul 


W  IIWU 


•  ••• 

J  Jf*t 
4  «««C 

—  mm  m  o 

O  0  00  7 
7  7  77  • 


5JJ  MIC  ••  M  C  o  o  o  • 

#  •  •  •  *o 

J\  JVK  O  QC  ft  -  ■  ..  ,  » 

5  r?  ^  7  22  2  ~  ~  ~  ^  •  •  •  •« 


•l  u  iluot  ft  ft  ft  o 

•£■  •>*  ^  •*«»««»- 
■  ►  >KIV  O  0  0(1  % 
—  >  O  ••••!) 


»-  (IOUU 

tt  •  •  •  «M 

m  —  nm  •  o 


•ft  OUftOU?  *-  K>  KK  5  % 

»U  Z  ^  u  ZK  •  J  J  J  j  *  %  ft  ft  UUUU  • 

ft  w»-  ^  wh  ••  •  •««  7  M  O  <  «  ««  ft  «• 

*►"  ooh-ooo  **N«niO  o  (i  o  k.^^h.o  «£ 


*s 

►  Al  AIX 
•  —  7 


l  -JfcJ  -4 tall* 

'lyoiio  v 

i  e  ♦  o  •  s 

•  •  •  03\33\D3\DD\  »  •  •  *2«  \  c"  • 

1*44044  04  404lL0ft  (HM  o  •  C  s  O  t  \ 

A^^^huu^^uwft^  ».oj  JJ4U&7S  u  "uu”gu"i-  "ft  Ji-ftiSli 

NWWMklODAwOOMCOO  •  •  •  O  •  O  •>  ■  Mi  ■  o  «f)  •  U  UUU  O  wft  Wft-ft» 

£z  X  i  x  so  m  c  z  z  z  —m  m  **  —  m  m  ««««<«  x  m  +  m  m  *. 

J<»  J«Z  *  I  ao  OOUJJftJJftJJftJJftkkkk  Iftl  ftS  s 


•  •0«w«C*  4444  1ft  Ift  t  ft  |ft  2 

•  444  44U  OU  <  «J  «U  •••*• 

£♦*£•«  OOOOft-O  £>*  1*0  »~0  2**2  J 

K  ^  t  ^  J  »  tu  •  uj  *  4i  «  o  o  o  o  •  •. 

•®  •  ••  *9  •  •  •  »x  #ftf  «ftt  Oftto  ft  ••  •  «ov 

••  *0(4  *07  4«4  4KOOK  O  OOft  O  «  Z  % 

'■**  <  0  0004  *0  4  *0  4  *0  4  *  O  •  •  •  •<  I 

III  u  ftIMU  ft  •  •••  •  *414  «l*l  ♦  •449000  •  ^ 

N04<wtt«HW04we«  «  I.  •  •  e 

v  ►  ft  wf>IV  ft  •  •  •  •  ftWK  WK  VW  ft  ft  K  ft«< 

x  >4  X  \4  c  OOOOZtwZfvx  <WZ2VJJJ.|C  ft 

*  ••  •  03\3D\D3\33\  •  •  •  «2  o 


V  rnm 
•  X  ft 

* 

4« 

4  ►*  *  2 
441  *0 
«X\M 
«  ► 

4  %(l  ft* 
4^*410 

4  m  *o 

4  *•  4 

4  ft* 

4  0  7  2 
«MO  M 
4  *K 
4  ftuy 
4  •-»  1*1 1- 
«*>< 
4  2  J 


4  ft! 

4  0  2 
MO 
4  *2 

v  c  * 
NO*  \ 


C  X  X  X  X 
2  4-  ►  J  J 

**  ■  ■  •  >  J  J 

m  04  o  O  O  O 


*  O  0C22I  ft*  M  «4  44444Z4 

J4  f  J4Z  40  tft 

j**jwwv^t*ww  vj  joj  jy  j  ju  j  jg  ww 


MftHUftHft f Skhmk ug^gu4U 


fOOOO*****.* 
Mft.  UUltHMHMft.  Z 


X  K  X  2 
X  ft*  7  2 

c  tc  o 

It.  K  It-  4 


9  CO 

—  Mft 

?SS2££22£22£  «oo  44 0  2  0  90 oooee ooo o  c  o  c  ooeeoeo 
e  e  •  44oo  o  o  eo  c  ccftiwwoir  o  oo  o  c  oeoeeeecoeooeoooo 
9C00040  e  44oe  Cftf-  ••  •*•  ftftfr  ooooe  ©  e  o  ooo**  mm  «  00  oo  oe 

4  M J.  4  N  *•-•  *.  *.  V  *  *.  *,  *  *  *.  *.  *.  *  «»  *.  *  •.  *.  *  *.  •,  «  *.  •»  *  N  *.  *  I.  t.  ft  •.  %  • 

Ji£SS2  5?S5f  Si?S5  55i?p,5*r<*  7b**7®o©©  9C**rvn«  *r 

r|T04442<  444444444  4  4  4  ^f*  K  K  4  ft  ft  fk  9  c  C4  4  B  C  <  C  C  2 


-  A>  •**  9  t r  «  9>  *  9  C  ft-  AK* 

SftftftftftftX  t  zAC  4C444490  070799  9  C  90  C  4  0  6  040  **ft  **•» 

ftft*««Mftft«ftwM^p.Mp,p,a.«,p>  ••  ••  •**  -  ft**  *>  ••  A  M  A»  A  A  MMMMMWMMV 


I 


It  It  ^  ^ 


A 

u 


a 

e 


2 

►- 

% 


a 

ft* 

x 


a 

0 


0  «  ^ 

ft*  X  Z 

M  •  • 

P*  «M 


X 

tol 


a 

(3 


a 

e 


a. 

e 


* 

w  ^ 

*■- 

5- 

to- 

H 

At 

2  — 

A) 

toi 

tot 

X 

•4 

to- 

►- 

to 

tt 

ft 

ft 

ft  p*  tol 

tol 

• 

• 

toi 

% 

-15  0 

« 

ft# 

-* 

ft 

-# 

IM  ft  w 

w 

to* 

to* 

% 

ft*  ft  »<• 

*4 

X 

z 

#- 

#- 

#-  % 

-  —  * 

z 

ft 

* 

ft  ft 

ft  >  » 

to 

to* 

w* 

•n 

•to 

X 

0 

—  to- 

to 

to- 

to- 

ft 

ft 

%# 

#-  ft 

ft*  to- 

•- 

tol 

tol 

*■*  W 

*« 

#*  ftto 

•* 

-1  ft 

*  «>  to 

tol 

to 

ft 

W 

►*  ft# 

•to  V* 

ft*  X 

4  •- 

-2  ft 

to 

• 

• 

ft  to 

x  or 

ft  or 

X  ft 

2  tot 

15  to 

to 

ft 

% 

toi  « 

to  * 

toi  « 

A  — 

a  u 

ft  ft  to*  #■» 

#- 

#- 

-#  #— 

X  #ft 

W 

x  ft 

o  z 

X  •  2  — 

to 

ftl 

to 

to  *n 

#ft*to 

ft  tol 

#-  ft 

2  to* 

»  ft*  to*  (5 

15 

15 

19 

#- 

*-  «ft# 

ft-  ft# 

*  w 

to-  — 

o 

ft  to  w 

w 

«0* 

w 

to' 

•n  #*  « 

x  m 

#ft  0 

X  ft' 

m  »* 

**  J  to* 

to* 

to* 

to* 

• 

W  KI  tol 

—  ft) 

—  to) 

•M  U 

•  15  ft  **  Z 

2 

X 

X 

♦ 

C  #*ft*  ♦ 

ft#  ♦ 

w  ♦ 

to-  — 

X  2 

2  2  ft  •  • 

•  ft 

ft*  — 1 

•n  » 

#- 

tol  ►-•)#- 

ft#  — 

0  #- 

ft#  • 

ft  ft-* 

•*  4  15  ft  ft* 

ft  to 

ftXI 

-  to 

ft- 

«  «  tol  #* 

2  *to 

tot  ito 

2  ►* 

Z  ft-*  15 

toll* 

mc» 

«  15 

ft- 

#• 

#-  *S  ♦  ft 

to  % 

♦  » 

ft  ft 

tol  to. 

to  ft  «« 

—  #- 

#*  #- 

ft-  ft- 

X 

X 

*n  x  ft*  #-m 

»#  »to 

—  tol 

to-  #ft 

C  •  M  ft  C  O 

o  c 

©  o 

o  © 

ft* 

«# 

»  o  -#  iM  *# 

ft  ft*  to* 

ito  -# 

X  ft* 

ft 

ft  X  ft  O  ft* 

c  «* 

o  -« 

O  ft* 

to* 

to- 

1*1  lO  5)  ft  to 

•  ft  X 

«  ft 

*  ft# 

-1  tol 

rw  to  o 

ft*  •* 

ft*  ft* 

ft*  ft* 

z 

2 

ft#  o  tol  »to  « 

•  «  to 

—  « 

•  « 

ft  -1  «-*  ■  to-  ft  • 

•  ft 

ft  ft 

ft  ft 

to* 

to  tol  »  w  #- 

©  #>  X 

ft##- 

o  tei 

U  15 

ftWUCto 

to  to 

to  • 

to  to 

2 

• 

to 

ft  to  #-  or  — 

*  **  #ft  to  in 

ft  « 

ft*  Z 

*  I  to  oo 

©  o 

o  o  «n  o  o  o 

O 

«* 

© 

#ftirt  ft*  «  ft# 

p*^ft 

•  -# 

— 

A-  -to 

—  to-  to  ft*  ft* 

ft-  ft* 

**  ** 

*4  K 

ft* 

X 

to 

to 

•n  —  w#-  m 

«ft  I 

—  0 

■  *# 

« 

nt  ft#  w 

ft#  ft* 

ft*  w  U  w  w  V) 

to- 

ft# 

to  toi  to  ft*  to  «  —to) 

-*  tol  ft#  to!  tol 

•to  — 

W  «  «  0.  v  to)  to* 

tol  tol 

U 

« 

to)  x 

X  to  -*  ki  ft  1 

to  •  2 

ft#  • 

«*  e 

> 

•  toto  to 

5-  to¬ 

to-  to- 

to  »-*-  to 

2 

© 

to 

ft-  tol  tin  ft«#* 

#«*M*  C  ^ 

to  to) 

0  U  0H  M  H 

Vto 

ft*  to* 

to 

to*  M 

to*  to*  « 

D 

• 

toi  to 

C  m  ft2  tol  — 

•  •to  « 

tol  •* 

•  ft 

to) 

Asortoa  stottakftsto 

to 

ft#  ^ 

e 

2 

«ti 

ft  2  X—  D  •  ft 

O  ft  ft 

•  • 

©  ft 

ft*  • 

«  ffzi iur iuii 

toi  2  2  tol 

X  X 

© 

•#*  X  x  O  ft*  ■  to  #-  to# 

.  •  toi  • 

#-  •» 

•  O 

X  w 

O  to  —  — 

4 

#-  #- 

4 

ft»  ft- 

to 

—  —  3 

Or 

•  ft* 

%# 

► 

M  X  NU  to)  ft# 

fti  ft#  • 

'•to  %# 

01  — 

to  to- 

It  rarvft  Kitft 

•  o  to  «  2 

toi 

X«9)-Z 

D 

2 

2  m  *.  m  ft  » 

•x  ft  e 

ft* 

•  2  ft# 

to 

toft  Z  OO 

«  o  c  « 

DOW 

c  o  to  a. 

(5 

a 

2 

to 

© 

O 

C  2  ©  tol  — 

©  •  a  c 

»to  w 

O  *0| 

•  T 

to  to*  2  2  to 

2  ?  to 

22 

k  z  z  )-r 

to  _Jto 

to* 

3 

or 

Z 

Zteift— -Jwt**-*-  — 

ft#  1 

ttoto 

ft  (K 

<5  0  ft*  w 

tol 

W  W 

tol 

ft#  ft* 

u 

wwy>-OI-<  m 

*  >-  O  C 

X 

Hi*  -j  to  a 

•  ^  » 

ft  0 

■  W  ft* 

«to  c 

»  to  to 

a 

to  to 

to 

to  to 

ft 

to  to  to  tol 

2  2 

(41  — 

toft 

tol  J 

P  D  C  Oft*  to  O  to  ■  *# 

-to  to) 

■  —•toto.  ft 

ft-  to 

ft  ft  ft  to*  to* 

ft 

to*  pi 

« 

to*  to* 

to 

to*  to*  to  to 

tol  to* 

to  x 

X 

to  n 

j  m 

u 

OOOCOt-ft-to—  to 

►  to- 

ft  to*  X 

to 

♦  ♦  ♦ 

♦ 

♦ 

♦ 

♦ 

oco900occoooeocoooec^oooooccoooecocoooccooo 
occcocooocecococccoccirtoeeccoe  o  o  c  c  c  &  e  oc  c  eoec 
ooeftoootrocinooiro  onto©— nj'vmft'roooooococeooec  s  o  c  o 

.  «  •.  I.  t  «,  «  «.  t  *.  t  I.  I,  t  t  ^  t  •.  t  ».  t  *.  «.  I.  t  •»  t  t  «.  «,  t  I.  ^  t  «.  «.  t  *.  t  t  I 

«  A-  «  t  O'  ©——  Alltofnftft»r.eto>toA*0  0  a  0  •  A  C  —  fV  m*H>«KCro#»Mf«»9l^3h. 


oo  e  c oooooc  cooe 


»t«,000  10  FORMATf*  FOR  ROTATION#  », J( *fc ( », II , •)  ■  *,6) 

919.000  ♦  »  ACTUAL  ANGLES  ARE*,/ 


C  oCOOCCwCOOC  sooooccccoooc 
o  c  e  ee  e  e  ec  e  e  c  c  oe  eec  o  c  o  c  e  c  c 
ir  4  odouwceoeoooeooocc  o  o  o  o  o 

I.  I.  i  A  L  I.  I.  A  i  L  I.  >■  L  A  t,  A.  L  *>  A  I.  I.  L  L  >  I 


ELECTRONICS  TECHNOLOGY  AND  DEVICES  LABORATORY 
MANDATORY  CONTRACT  DISTRIBUTION  LIST 


101  Defense  Technical  Information  Center 
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Cameron  Station  (Bldg  5) 

012  Alexandria,  VA  22314 

203  GIDEP  Engineering  &  Support  Dept 
TE  Section 
PO  Box  398 

001  NORCO,  CA  91760 

205  Director 

Naval  Research  Laboratory 
ATTN:  CODE  2627 
001  Washington,  DC  20375 

301  Rome  Air  Development  Center 
ATTN:  Documents  Library  (TILD) 

001  Griffiss  AFB,  NY  13441 

437  Deputy  for  Science  &  Technology 
Office,  Asst  Sec  Army  (R&D) 

001  Washington,  DC  20310 

438  HQDA  (DAMA-ARZ-D/Dr.  F.D.  Verderame) 
001  Washington,  DC  20310 
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US  Army  Materiel  Systems  Analysis  Actv 
ATTN:  DRXSY-MP 
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001  Alexandria,  VA  22333 
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001  Adelphi,  MD  20783 

612  Cdr,  ERADCOM 
ATTN:  DRDEL-CT 
2800  Powder  Mill  Road 
001  Adelphi,  MD  20783 

680  Commander 

US  Army  Electronics  R&D  Command 
000  Fort  Monmouth,  NJ  07703 
1  DELET-MQ 
1  DELEW-D 
1  DELET-DD 

1  DELSD-L  (Tech  Library) 

2  DELSD-L-S  (STINFO) 

34  Originating  Office 

1  DELET-MF 

681  Commander 

US  Army  Communications  R&D  Command 
ATTN:  USMC-LNO 
001  Fort  Monmouth,  NJ  07703 


ELECTRONICS  TECHNOLOGY  AND  DEVICES  LABORATORY 
SUPPLEMENTAL  CONTRACT  DISTRIBUTION  LIST 
(ELECTIVE) 


103  Code  R123,  Tech  Library 
DCA  Defense  Comm  Engrg  Ctr 
1800  Wiehle  Ave 

001  Reston,  VA  22090 

104  Defense  Communications  Agency 
Technical  Library  Center 

Code  205  (P.  A.  Tolovi) 

001  Washington,  DC  20305 

206  Commander 

Naval  Electronics  Laboratory  Center 
ATTN:  Library 
001  San  Diego,  CA  92152 

207  Cdr,  Naval  Surface  Weapons  Center 

White  Oak  Laboratory 

ATTN:  Library  Code  WX-21 
001  Silver  Spring,  MD  20910 

314  Hq,  Air  Force  Systems  Command 

ATTN: DLCA 
Andrews  Air  Force  Base 
001  Washington,  DC  20331 

403  Cdr,  MICOM 

Redstone  Scientific  Info  Center 
ATTN:  Chief,  Document  Section 
001  Redstone  Arsenal,  AL  35809 

406  Commandant 

US  Army  Aviation  Center 
ATTN:  ATZQ-D-MA 
001  Fort  Rucker,  AL  36362 

407  Director.  Ballistic  Missile  Defense 
Advanced  Technology  Center 
ATTN:  ATC-R,  PO  Box  1500 

001  Huntsville,  AL  35807 

418  Commander 

HQ,  Fort  Huachuca 
ATTN:  Technical  Reference  Div 
001  Fort  Huachuca,  A Z  85613 


475  Cdr,  Harry  Diamond  Laboratories 
ATTN:  Library 
2800  Powder  Mill  Road 
001  Adelphi,  MD  20783 

477  Director 

US  Army  Ballistic  Research  Labs 
ATTN:  DRXBR-LB 

001  Aberdeen  Proving  Ground,  MD  21005 

*481  Harry  Diamond  Laboratories 

ATTN:  DELHD-RCB  (Dr.  J.  Nemarich) 

2800  Powder  Mill  road 
001  Adelphi,  MD  20783 

482  Director 

US  Army  Materiel  Systems  Analysis  Actv 
ATTN:  DRXSY-T,  MP  (In  Turn) 

001  Aberdeen  Proving  Ground,  MD  21005 

507  Cdr,  AVRADCOM 
ATTN:  ORSAV-E 
PO  Box  209 

001  St.  Louis,  MO  63166 

511  Commander,  Picatinny  Arsenal 
ATTN:  SARPA-FR-5,  -ND-A-4, 

-TS-S  (In  Turn) 

001  Dover,  NJ  07801 

515  Project  Manager,  REMBASS 
ATTN:  DRCPM-RBS 
001  Fort  Monmouth,  NJ  07703 

517  Commander 

US  Army  Satellite  Communications  Agcy 
ATTN:  DRCPM-SC-3 
001  Fort  Monmouth,  NJ  07703 

518  TRI-TAC  Office 
ATTN:  TT-SE 

001  Fort  Monmouth,  NJ  07703 


For  Millimeter  &  Microwave  Devices  Only 


ELECTRONICS  TECHNOLOGY  AND  DEVICES  LASORATORY 
SUPPLEMENTAL  CONTRACT  DISTRIBUTION  LIST  (CONT) 
(ELECTIVE) 


519  Cdr,  US  Army  Avionics  Lab 
AVRADCOM 

ATTN:  DAVAA-D 

001  Fort  Monmouth,  NJ  07703 

520  Project  Manager,  FIREFINDER 
ATTN:  DRCPM-FF 

001  Fort  Monmouth,  NJ  07703 

521  Commander 

Project  Manager,  SOTAS 
ATTN:  DRCPM-STA 
001  Fort  Monmouth,  NJ  07703 

531  Cdr,  US  Army  Research  Office 
ATTN:  DRXRO-PH  (Dr.  lontz) 

DRXRO-IP  (In  Turn) 

PO  Box  12211 

001  Research  Triangle  Park,  NC  27709 

556  HQ,  TCATA 

Technical  Information  Center 
ATTN:  Mrs.  Ruth  Reynolds 
001  Fort  Hood,  TX  76544 

568  Commander 

US  Army  Mobility  Eqp  Res  &  Dev  Cmd 
ATTN:  DRDME-R 
001  Fort  Belvoir,  VA  22060 

604  Chief 

Ofc  of  Missile  Electronic  Warfare 
Electronic  Warfare  Lab,  ERADCOM 
001  White  Sands  Missile  Range,  NM  88002 

606  Chief 

Intel  Materiel  Dev  &  Support  Ofc 
Electronic  Warfare  Lab,  ERADCOM 
001  Fort  Meade.  MD  20755 


608  Commander 
ARRADCOM 
DRDAR-TSB-S 

001  Aberdeen  Proving  Ground,  MD  21005 

614  Cdr.  ERADCOM 

ATTN:  DRDEL-LL,  -SB.  -AP  (In  Turn) 
2800  Powder  Mill  Road 

001  Adelphi,  MD  27083 

617  Cdr,  ERADCOM 
ATTN:  DRDEL-AQ 
2800  Powder  Mill  Road 

001  Adelphi,  MD  20783 

619  Cdr,  ERADCOM 

ATTN:  DRDEL-PA,  -ILS.  -ED  (In  Turn) 
2800  Powder  Mill  Road 

001  Adelphi,  MD  20783 

701  MTI  —  Lincoln  Laboratory 
ATTN:  Library  (RM  A-082) 

PO  Box  73 

002  Lexington,  MA  02173 

703  NASA  Scientific  &  Tech  Info  Facility 
Baltimore/Washington  (ntl  Airport 

001  PO  Box  8757,  MD  21240 

704  National  Bureau  of  Standards 
Bldg  225,  RM  A-331 

ATTN:  Mr.  Leedy 

001  Washington,  DC  20231 

707  TACTEC 

Batelle  Memorial  Institute 
505  King  Avenue 

001  Columbus,  OH  43201 
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Coordinated  Science  Laboratory 
University  of  Illinois 
Urbana,  Illinois  61801 

ATTN:  Dr.  Bill  J.  Hunsinger  (1) 

Dr.  J.S.  Bryant 
OCRD 

ATTN:  DARD-ARP 

Washington.  DC  20310  (1) 

Dr.  R.  LaRosa 

Hazeltine  Corporation 

Greenlawn,  New  York  11740  (1) 

General  Electric  Co. 

Electronics  Lab 
Electronics  Park 
Syracuse,  NY  13201 

ATTN:  Mr.  S.  Wanuga  (1) 

Air  Force  Cambridge  Labs 
ATTN.  CRDR  (Dr.  P.  Carr  &  Dr.  A.J. 
Slobodnik) 

Bedford.  MA  01730  (2) 

Mr.  R.  Weglein 

Hughes  Research  Laboratories 

3011  Malibu  Canyon  Road 

Malibu,  California  90265  (1) 

Mr.  H.  Bush  CORC 
RADC 

Griffiss  Air  Force  Base 

New  York  13440  (1) 

Dr.  Tom  Bristol 
Hughes  Aircraft  Company 
Ground  Systems  Group 
Bldg  600/ MS  D235 
1901  W.  Malvern 

Fullerton,  CA  92634  (2) 


Anderson  Laboratories,  Inc. 

1280  Blue  Hills  Ave 
ATTN:  Dr.  A.A.  Comparini 
Bloomfield,  Conn.  06002 

Mr.  Henry  Friedman 

RADC/OCTE 

Griffiss  AFB,  NY  13440 

Autonetics,  Division  of  North 
American  Rockwell 
P.0.  Box  4173 
3370  Miraloma  Avenue 
Anaheim,  CA  92802 
ATTN:  Dr.  G.R.  Pulliam 

General  Dynamics,  Electronics  Division 
P.0.  Box  81127 
San  Diego,  CA  92138 
ATTN:  Mr.  R.  Badewitz 

Texas  Instruments,  Inc. 

P.0.  Box  5936 

13500  N.  Central  Expressway 
Dallas,  Texas  75222 
ATTN:  Dr.  L.T.  Clairborne 

Raytheon  Company 
Research  Division 
28  Seyon  Street 
Waltham.  Massachusetts  02154 
ATTN:  Dr.  M.B.  Schulz 

Sperry  Rand  Research  Center 
100  North  Road 

Sudbury,  Massachusetts  01776 
ATTN:  Or.  H.  Van  Oe  Vaart 

Microwave  Laboratory 

W.W.  Hansen  Laboratories  of  Physics 

Stanford  University 

Stanford,  CA  94305 

ATTN:  Dr.  H.J.  Shaw 


Commander,  AFAL 

ATTN:  Mr.  W.J.  Edwards,  TEA 

Wright-Patterson  AFB,  Ohio  45433  (1) 
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Polytechnic  Institute  of  Brooklyn 

Route  No.  110 

Farmingdale,  NY  11735 

ATTN:  Dr.  A.A.  Oliner  (1) 

Westinghouse  Electric  Corp. 

Research  &  Development  Center 
Beulah  Road 
Pittsburgh.  PA  15235 

ATTN:  Dr.  J.  DeKlerk  (1) 

Stanford  Research  Institute 

Menlo  Park,  CA  94025 

ATTN:  Dr.  A.  Bahr  (1) 

International  Business  Machines  Corp. 
Research  Division 
P.O.  Box  218 

Yorktown  Heights,  NY  10598 

ATTN:  Dr.  F.  Bill  (1) 

TRW 

Defense  and  Space  Sys  Group 

One  Space  Park 

Redondo  Beach,  CA  90278 

ATTN:  Dr.  R.S.  Kagiwada  (1) 

Tektronix  Inc. 

P.O.  Box  500 
Beaverton,  OR  97077 

ATTN:  Dr.  R.  Li  (1) 

Dr.  Fred  S.  Hickernell 

Integrated  Circuit  Facility 

Motorola  Government  Electronics  Division 

8201  East  McDowell  Road 

Scottsdale.  AZ  85257  (1) 

Prof.  H.F.  Tiersten 
Jonsson  Engineering  Center 
Rensselaer  Polytechnic  Institute 
Troy.  NY  12181  (1) 

McGill  University 
ATTN:  G.W.  Farnell 
Montreal  110,  Canada 


Advanced  Technology  Center.  Inc. 
Subsidiary  of  LTV  Aerospace  Corp. 

P.O.  Box  6144 
Dallas,  Texas  75222 

ATTN:  Mr.  A.E.  Sobey  (1) 

United  Aircraft  Research  Labs 

ATTN:  Dr.  Thomas  W.  Grudkowski 

East  Hartford,  Conn.  06108  (1) 

Science  Center 

Rockwell  International 

Thousand  Oaks,  CA  91360 

ATTN.  Dr.  T.C.  Lim  (1) 

University  of  Southern  CA 
Electronic  Science  Lab 
School  of  Engineering 
University  Park,  Los  Angeles 
California  900 

ATTN:  Or.  K.  Lakin,  SSC  303  (1) 

SAWTEK,  Inc. 

P.O.  Box  7756 

2451  Shader  Road 

Orlando.  Florida  32854 

ATTN:  Mr.  S.  Miller  (1) 

Prof.  P.C.Y.  Lee 
School  of  Engineering 
Princeton  University 

Princeton.  NJ  08540  (1) 

Mr.  John  A.  Kusters 

Hewlett-Packard 

5301  Stevens  Creek  Boulevard 

Santa  Clara.  CA  95050  (1) 

Dr.  Tom  Young 
Sandia  Laboratories 
P.O.  Box  5800 

Albuquerque,  NM  87185  (1) 

Dr.  William  J.  Tanski 
Sperry  Research  Center 
100  North  Road 
Sudbury,  MA  01776 


(D 


(D 
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Dr  B.A.  Auld 

Dr.  Robert  L.  Rosenberg 

E  L.  Ginzton  Laboratory 

Bell  Laboratories 

Stanford  University 

Crawfords  Corner  Road 

Stanford.  CA  94305 

(D 

Holmdel.  NJ  07733 

(D 

Mr.  Marvin  E.  Frerking 

Dr.  B.K.  Sinha 

MS  137-138 

Schlumberger-Doll  Research  Center 

Collins  Radio  Company 

P.O.  Box  307 

855  35th  Street.  NE 

Cedar  Rapids.  IA  52406 

(D 

Ridgefield.  CT  06877 

Dr.  Robert  S.  Wagers 

(D 

Dr.  William  R.  Shreve 

Texas  Instruments,  Inc. 

HP  Laboratories 

13500  N.  Central  Expwy. 

1501  Page  Mill  Road 

P.O  Box  225936.  MS  134 

Palo  Alto.  CA  94304 

(1) 

Dallas.  TX  75265 

(D 

Dr.  Thomas  M.  Reeder 

Dr.  Richard  C.  Williamson 

MS  50-362 

Lincoln  Laboratory 

Tektronix.  Inc. 

P.O.  Box  73 

P.O.  Box  500 

Beaverton,  OR  97077 

(D 

Lexington,  MA  02173 

0) 

SUPPLEMENT  TO  DISTRIBUTION  LIST 

D.  Chrissotimos.  Code  763 
National  Aeronautics  and  Space 
Administration 
Goddard  Space  Flight  Center 
Greenbelt.  MD  20771  (1) 

Naval  Research  Laboratories 
Code  5237 

Washington.  DC  20375 
ATTN:  Dr.  D.  Webb  (1) 

HQ  ESD  (DRI) 

L.G  Hanscom  AFB 

Bedford.  MA  01731  (1) 

Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-RE  (Mr.  Pittman) 

Redstone  Arsenal,  AL  35809 


Army  Materials  and  Mechanics  Research 
Center  (AMMRC) 

Watertown,  MA  02172 

ATTN:  DMXMR-EO  (1) 

Commander,  Picatinny  Arsenal 
ATTN:  SARPA-FR-S 
Bldg.  350 

Dover,  NJ  07801  (2) 

A.  Kahan 
RADC/ESE 
Hanscom  AFB 

Bedford,  MA  01731  (1) 


(D 


